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The  stochastic  duel  was  first  formulated  and  analyzed  by  G. 
Trevor  Williams  and  the  author,  in  a paper  presented  on  11  May  1966 
at  the  k-th  Annual  ORSA  meeting  in  Washington,  D.C.  Since  then, 
about  fifty  research  papers  have  been  written  by  authors  scattered 
in  many  parts  of  the  world.  This  fairly  substantial  body  of  results 
has  proven  of  value  to  many  military  OR  analysts,  both  in  and  out  of 
the  services.  However,  much  of  it  is  inaccessible,  or  mathematically 
obtuse.  Consequently,  we  have  set  out  to  remedy  this  situation,  with 
the  support  of  the  U.  S.  Amy  Research  Office. 

This  report  is  approximately  one-half  of  the  total  output  of 
the  project.  The  material  prepared  so  far  is  sufficiently  self- 
contained  as  to  warrant  issuance  at  this  time.  The  remainder  of  the 
report  will  be  forthcoming  within  the  next  year,  at  which  time,  the 
entire  manuscript  will  be  updated  and  brought  together  as  one  entity. 

The  present  work  consists,  primarily,  of  two  parts.  In  Part 
I an  exposition  of  the  two  principal  methods  of  deriving  results  is 
given.  These  are  the  mixture  technique  and  the  semi-Markov  terminat- 
ing renewal  process  technique.  For  the  first  time,  to  the  best  of 
the  author’s  knowledge,  it  is  carefully  and  explicitly  shown  that, 
in  fact,  these  are  the  techniques  being  used  and  precisely  how  they 
are  being  used;  and,  finally,  how  they  are  related.  We  hope  to 
clear  up  much  of  the  obscurity  and  mystery  surrounding  the  utility 
of  the  techniques  and  also  to  i:  Ucate  when  there  is  (or  is  not)  an 
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advantage  to  using  one  or  the  other.  This  is  done,  primarily,  by  carry- 
ing an  example  through  each  technique. 


I 


I 


/ 

\ 


Bart  II  is  a comprehensive,  exhaustive  and  fully  annotated  bib- 
liography of  all  research  papers  on  the  topic,  known  to  the  author. 

It  has  been  formatted  to  uniformly  display,  in  detail,  what  has  been 
accomplished  in  each  paper. 

The  additional  work,  to  be  done  later,  will  include: 

(1)  an  historical  and  expository  introduction; 

(2  ) possibly,  some  further  exposition  on  techniques,  and, 
most  importantly, 

(3)  a comprehensive  compendium  of  all  known  results  reduced 
to  a common  notation  and  organized  to  facilitate  the 
location  of  any  desired  result. 

It  is  believed  that  this  project  will 

(1)  make  all  results  easily  and  conveniently  available  to  the 
analyst,  and 

(2 ) will  aid  the  research  worker  in  identifying  what  has  been 
done  (to  prevent  reinventing  the  wheel)  and  what  has  not 
been  done,  so  that  he  may  direct  his  efforts  in  a produc- 
tive direction. 

The  author  would  be  grateful  to  learn  of  any  inadvert  ant  mis- 
takes, omissions,  or  other  errors  which  may  have  occurred. 
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PART  I 

AN  EXPOSITION  OH  TECHNIQUES 

1.  INTRODUCTION 

In  the  fundamental  duel  (FD),  two  contestants , A and  B,  fire 
at  each  other  at  certain  interval*  and  either  hit  ox  alee  on  each  round 
fired.  The  duel  teminatee  when  either > or  both,  are  hit.  The  hit 
probabilities  for  each  are  constant  from  round  to  round  and,  in  general, 
are  different  for  each.  The  time  between  rounds  nay  be  a continuous  random 
variable  or  may  be  a constant,  and,  in  general,  is  different  for  each.  The 
contestants  both  start  at  time  sero  with  unloaded  weapons  and  fire  their 
first  rounds  sometime  later.  They  both  have  unlimited  ammunition  and  un- 
limited time. 

The  analysis  begins  by  first  considering  a marksman  firing  at  a 
passive  target,  under  the  same  conditions  given  above,  until  he  hits  it. 

This  is  called  the  fundamental  marksman  problem  (FN).  From  this  we  may 
solve  the  duel  problem  by  considering  each  marksman  to  be  firing  independ- 
ently of  the  other.  The  first  to  hit  his  passive  target  wins.  This  is 
entirely  equivalent  to  the  fundamental  duel,  as  the  model  described  above 
in  no  way  links  the  actions  of  one  duelist  to  those  of  the  other. 

Consequently,  let  us  first  consider  the  case  of  the  marksman  versus 
a passive  target,  and  further,  let  vs  confine  ourselves  here  to  the  situa- 
tion where  his  Interfiring  time  is  a continuous  random  variable.  It  is 
clear  that  he  might  hit  on  his  first  round  fired,  or  possibly  on  the 
second,  or,  in  fact,  on  any  round,  providing  he  has  failed  to  hit  on  all 
preceding  rounds.  That  is,  if  his  hit  probability  is  p,  (q  ■ 1-p),  then 
the  probability  that  his  first  (and  fatal)  hit  is  an  the  nth  rexind  and  is 
p q If  his  interfiring  time  has  % probability  density  function  (pdf) 
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given  by  f(t),  then  f(t)  la  also  hia  p df  to  a hit  on  the  flrat.  round. 
If  be  hlta  on  the  ascend  round  after  failing  an  the  flrat,  hla  pdf  to  a 
hit  la  the  convolution  of  f(t)  with  itself,  denoted  by  f*f(t).  Thla  la 
because  if  he  has  fired  twice,  he  haa  made  two  selections  (at  random)  from 
f(t)  and  added  then  together  to  determine  hla  tine  to  fire  the  eecond 
round.  Continuing,  hla  pdf  of  tine  to  fire  the  round  la  given 

by  f*f*f*  ...  *f(t)  • f^t),  i.e.,  n convolutions  of  f(t)  with 

itself. 

The  stochastic  process  described  above  la  called  a mixture  and  may  be 
described  as  follows.  Let  T be  the  continuous  random  variable,  tine  to  a 
hit,  and  let  X^  1 ■ 1,2,3,...,  be  the  continuous  random  variables,  tine 
between  firing  epochs,  where  the  are  known,  independent  and  identically 
distributed  (lid)  as  X with  pdf  f(t).  Then 

T - Xx 

-Xl+  *2 

e 

" *1  + £,  + X5  + •••  + xn 

e 

e 

where  all  the  rows  above  are  the  mutually  exclusive  end  exhaustive  ways  of 
obtaining  a hit.  Now,  let  h(t)  be  the  pdf  of  T and  f(t)  be  the 
of  each  X^,  then  by  the  basic  property  of  mixtures: 


with  probability 


P 

pq 


pq 


> (l) 


n-l 


pdf 


■ 


or 


P[t  < T < t ♦ dt) 


■ h(t)dt  » pf(t)dt  ♦ pqf2#(t)dt 
♦ pq2f^*(t)dt  ♦ •••  ♦ pqn-ljP* (t )dt  ♦ .. 


h(t) 


^ pqn_1 


n«l 


(2) 


This  srprsssion  can  be  greatly  simplified  by  converting  each  side  to  a 
characteristic  function  (Fourier  transform)  as  follows : 


u 

i 


*(u)  - [ eiut  h(t)dt  and  ?(u)  - f eiut  f(t)dt  . (3) 

J0  J0 

Using  the  definitions  (3)  in  Equation  (2),  and  the  convolution  property  of 
characteristic  functions  (cf): 


♦(u) 


I 

n-1 


n-1  \ 

q 9 (u) 


■ p<P(u)  [q<p(u)]n-1 


n-1 


_.L£&L_  , 

1 - q <p(u ) * 


00 


where  the  indicated  sum  is  valid  because  qf(u)  < 1,  for  all  u.  Equation 
(k)  can  now  be  inverted  to  give 


h(t)  " 2Ti  J *(u)du 


[m  e"1Ut  yCuldu 
2 n J „ 1 - q*(u} 


(5) 
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We  should  note,  la  passing,  that  the  result  in  Equation  (4)  can  be  written 
down  at  cnee  after  inspecting  Equation  (l).  This  cones  about  because  In 
(l),  the  sequence  of  probabilities  p , pq  , . . . , pq®"1  , . . . is  the  proba- 
bility ease  function  (pnf ) of  the  randan  variable  X,  the  round  number 
on  which  a hit  occurs.  It  is  geonstric  and  given  by 

Pjj(b)  ■ W » n * 1,2,..., 

■ 0 , elsewhere  , 

which  has  a geometric  transform 

m m 

Gj|(z)  ■ Pj,(n)*n  " P qn_1  *n 

n * - “ n«l 


n»l 


There  is  a theorem  (see,  Giffla,  1975,  Eqn.  4-28)  which  says  that  for  a 
mixture  of  lid  random  variables: 

•(u)  - GjjtfCu))  . (7) 

Therefore,  from  (5)  and  (6),  we  have  immediately  $(u)  • » as 

before . 

The  solution  to  the  duel  may  now  be  written  deem.  If  ?A  and  T£ 
are  the  random  variables,  times  for  A and  B to  hit  a passive  target, 
with  pdf's  hA(t)  and  hg(t),  respectively,  then  the  probability  that 
A wins  the  duel,  P[A],  is 
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P[A]  - P[Ta  < Tb) 


/ P[t  < Ta  < t + dt.  , Tb  > t]dt 

J hA(t)  }£(t)dt  , 


(8) 


where  the  third  line  is  Justified  by  the  Independence  of  TA  and  Tfi. 

This  nay  he  put  in  more  tractable  form  using  characteristic  func- 
tions by  a theorem  of  Barseval  which  states. 


J 


f^x)  f2(x)dx 


^(-u)  (u  )du  . 


Using  (9)  in  (8), 


(9) 


P[A] 


V-u) 


[•B(u)  - 1] 

— 


du 


(10) 


It  is  easily  shown  that  the  ®(-u)  function  has  no  poles  in  the  complex 
lower  half-plane  and  that  it  vanishes  on  a large  semicircle,  C,  in  the 
lower  half-plane  as  R -♦  ••  (see  Fig.  l). 

Thus,  if  we  go  to  complex  u,  by  analytic  continuation,  the  second 
term  in  the  integrand  is  xero  and  we  have 
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p[a3 


ie 

ie 


®A(-u)  *B(u)du 
u ““ 


(n) 


where  6 is  less  than  the  distance  to  the  nearest  pole  in  the  lower  half 
of  the  complex  u plane.  Since  a draw  is  impossible  with  continuous 
firing  times,  P(B]  - 1 - P[A]. 


• ^sC5Z  ^ 


Much  of  the  literature  uses  the  Laplace  transform  (LT)  of  the  time 
functions  involved,  hut  we  shall  always  use  the  characteristic  function  far 
two  reasons : 

(1)  the  cf  always  exists,  whereas  the  LT  may  not,  and 

(2 ) the  numerical  Integration  of  the  Fourier  transform  (cf  ^ 
has  received  extensive  attention  and  many  efficient 
algorithms  are  available  (see,  for  example,  Brigham,  197^ ) . 

Finally,  let  us  emphasize  the  efficiency  and  economy  of  using  charac 
teristic  functions  by  comparing  Equation  (l)  with  the  solution  to  the  duel 
without  their  use.  Hitting  (2)  with  appropriate  subscripts  for  A and  B 
into  (8),  we  have 


P[A] 


/0~  z fT<t)v  / ' i 


n=l 


t n»l 


and  from  (ll). 


P[A] 


JL  f"-1 

2Ti  J m m_ 


9„(u) 


*-ie  *A  JbV-  du 

ie  [1  " *a  [1  * *B  u 


SH  . ( 


We  see  that  in  (13)  we  have  one  integration  to  perform,  whereas  in  (12 ),  we 
have  two  integrations  of  two  infinite  sums,  each  term  of  which  Involves 
iterated  convolutions.  This  is  indeed  an  enormcus  simplification. 


2.  SOLUTION  METHODS 


In  this  section  we  shall  examine  several  solution  techniques  commonly 
used  in  the  literature.  For  this  purpose,  let  us  solve  the  FM  problem, 
slightly  modified  as  follows.  A marksman  fires  with  constant  hit  probabil- 
ity, p,  at  a passive  target.  His  target  is  destroyed  when  he  hits  it 
twice.  His  interfiring  times  are  continuous,  lid,  and  have  a general 
(unspecified)  pdf. 

A.  The  Mixture  and  Characteristic  Function  Technique 

Using  the  notation  given  above  and  proceeding  as  before, 


T SX1+X2 
- Wx3 
-V  WXU 

” W Xn 

e 

Mote  that,  in  general,,  the  nth  term  has  a coefficient  ( n ) because 

we  can  have  the  first  hit  anywhere  in  the  first  n-1  positions,  while  the 

th 

second  hit  must  be  on  the  n firing. 

Continuing  as  before, 

b(t)  - Y ( V’1  ) p2  q”-2  <**(0)  - 


with  probability 


V") 


2 


pqp  + qpp 

pqqp+  qpqp+  qqpp 


n-1 


2 n-2 


)p‘  q 


(1 


- I ( ;:i ) P2  ,n-2  , 

n^ 

and  talcing  the  cf  of  both  sides 

m 

•(«)  - <p2(u)P2  y < si 2 ><«*<«> )“'2 

•§ 

- [p?(u)l2  Y ( n*2*1  )U<PW]“  • (15) 

n=0 

Recognizing  the  stun  as  the  negative  binomial  series  with  parameter  2,  ve 
may  immediately  write 


*(u)  - p2  ?Z(u)(l  -q9(u)]‘Z 


L&fe).  . 

1 - q?(uj 


(16) 


Two  comments  are  in  order  here.  First,  the  pmf  for  the  round  number  (n) 
on  which  the  process  terminates  is  given  by  the  right-hand  column  of  (l4) 
and  is,  of  course,  a negative  binomial  with  parameter  2,  identical  to  the 
one  above.  Its  geometric  transform  is 


and  using  (7), 


» 


*(u) 


ZSfel. 

1 - q 9(u) 
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which  is  the  same  as  (l6).  Secondly,  we  ncte  that  this  result  is  obvious 
because  the  time  to  a second  hit  is  simply  the  time  to  a first  hit  plus  the 
independent  time  to  a second  hit,  both  of  which  are  identically  distributed 
as  the  time  to  a first  hit  whose  cf  is  given  by  (4),  and  therefore,  our 
result  follows  immediately.  Although  the  result  is  obvious,  it  is  still 
instructive  to  derive  it  by  several  different  techniques. 

B.  The  Semi-Markov  Process  Technique 

In  this  technique,  Markov  chain  theory  is  used  to  establish  the  pmf 
of  N,  the  round  number  on  which  the  contest  terminates . This  is,  of 
course,  the  right-hand  column  of  (l4)  and  one  might  wonder  why  it  would  be 
advantageous  to  use  a relatively  complicated  technique  such  as  this  to 
establish  what  appears  to  be  rather  obvious . The  answer  is  two- fold: 

(l)  to  demonstrate  the  technique,  and 

(2  ) when  the  state  space  is  more  complicated  than  in  our  example, 
the  derivation  can  become  very  difficult  as  there  will  be  so 
many  branches  in  the  outcome  tree  that  keeping  track  of  all 
the  paths  through  it  may  become  nearly  impossible.  The 
Markov  chain  does  this  automatically. 

The  first  step  is  to  establish  the  state  space,  which  in  our  case, 
has  three  states:  (l)  not  hit  (ft),  (2)  hit  but  not  killed  (hR),  i.e., 

one  hit  only,  and  finally,  (3)  hit  and  killed,  or  simply  killed  (K), 
i.e.,  a second  hit. 

We  notice  that  in  the  language  of  Markov  chains,  two  states  are 
transient  (HX , ft)  and  one  (k)  is  absorbing  (once  entered,  it  may 
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never  be  left).  Next,  the  initial  state  probability  vector  (i)  and  the 
transition  probability  matrix  on  the  state  space  (s)  are  determined.  In 


our  example  ve  have 


K HR  H 

(0,0,1) 


which  are  now  partitioned  as  follows : 


I = (0,  :o,  1)  - (0,  M) 


1 ! 0 


Oj  p 


1 0 


T P 


The  row  vector  M excludes  the  absorbing  state  K.  The  sub-matrix  P 
contains  all  the  transition  probabilities  for  remaining  in  non-absorbing 
states  and  the  column  vector  T contains  all  the  transition  probabilities 
for  going  from  a transient  state  (HK  or  S)  to  the  absorbing  state  K. 
From  elementary  Markov  theory  we  have  that  Pjj(n),  the  probability  of  a 
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kill  on  the  n'“  round  fired  is 


Vn)  ■ a?’1 1 • 


First,  determine  by  matrix  multiplication: 


*-c  :> 


a :■) 


• ° i 

3p42  ' 


( \ ° ) 

\ n-1  n J 


The  general  result  for  rf1  is  easily  proven  by  induction.  Continuing, 


M P0-1  - (0,1 


, , \ n-2  n< 

(n  - l)pq  q 


- ((n-l^.q^J 


-1  ) 


=1 


VB)  - E?'1 1 


- g ) 

- (n  - l)Mn"2  - x j m”"2  , 

which  is,  of  course,  the  general  term  of  the  pmf  of  H and  corresponds 
to  the  general  term  of  the  right-hand  column  of  (lU)  above  in  A.  Prom 
here,  the  development  follows  A above  precisely  and  will  not  be 
repeated.  The  reason  that  this  is  called  a semi- Markov  process  is  because 
the  time  between  events  is  a RV. 

Before  turning  to  the  next  topic,  we  should  note  that  the  state 
space  may  be  much  larger.  There  may  be  many  transient  states,  thus,  P, 

M and  T may  be  such  larger.  In  fact,  there  may  be  quite  a few  absorbing 
states.  For  example,  it  might  be  possible  to  be  killed  on  the  first  hit,  or 
not  killed  on  the  first  hit  and  killed  on  the  second  hit  (as  above).  The 
point  is,  that  this  technique  provides  a convenient  and  orderly  way  to 
expand  the  state  space. 

C.  Renewal  Theory  Integral  Equations  Technique 

It  will  be  helpful  in  looking  at  the  techniques  to  be  described 
in  what  follows,  to  establish  first  seme  notation  and  concepts  from 
renewal  theory. 

Let  N(t)  be  the  number  of  firings  up  to  time  t.  Clearly,  for 
every  fixed  t,  this  is  a discrete  RV.  The  situation  may  be  depicted 
graphically  by  a possible  realization  as  shown  in  Figure  2.  This  process 
is  called  a terminating  semi- Harkov  renewal  process.  Most  of  the  renewal 
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theory  literature  concerns  itself  with  nonterminating  processes,  but  we  may 

y 

still  profitably  consider  this  approach. 

In  the  ordinary  nonterminating  renewal  process,  one  important  result 
is  that  E[N(t)3  can  be  expressed  as  a simple  integral  equation.  Ir  our 
case,  however,  ve  are  more  interested  in  the  function  h(t)  which  can 
also,  in  a similar  manner,  be  derived  from  a series  of  integral  equations  in 
sequential  order.  In  order  to  do  this  we  define  three  functions,  as 
follows. 

First,  let  Tn  ■ Xa  be  the  RV,  time  to  nth  firing,  with 
pdf  ^(t),  then, 

h°(t)dt  * P[t  < T < t+  dt  , 0 hits  in  n trials] 
n n 

- P[t  < Tn  < t + dt  | 0 hits  in  n trials] 

• P[0  hits  in  n trials] 

■ f***(t)qn  dt,  n ■ 1,2, ...»  (l8) 

h^(t)dt  * Pit  < Tr  < t + dt  1 hit  in  n trials] 

■ P[t  < T < t+  dt  I 1 hit  in  n trials] 

n 

• P[l  hit  in  n trials] 

» f"*(t)(  J )pqn-1  dt,  n * 1,2 , . . . , (19) 

where  the  factor  ( ” ) is  necessary  because  the  hit  may  occur  on  any 
round  fired,  and 

» 
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h (t)  ■ P[t  < T < t+  dt , 1 hit  in  n-1  trials  and 
n n 

2nd  hit  on  the  trial] 

* P[t  < T < t + dt  | 1 hit  in  n-1  trials  and 
n 

2nd  hit  on  the  n**1  trial] 

• P[1  hit  in  n-1  trials  and  2nd  hit  on  the  nth  trial] 

- f**(t)  ( n'1  h2  q11"2  dt,  n-2,3,...,  . (20) 

Now,  notice  that  if  we  consider  the  sequence  n ■ 1,2,...,  and  consider 
the  sum  on  n that  (l8)  and  (19)  are  not  proper  density  functions,  but  that 
(20 ) is  (it  is  Just  the  mixture  defined  in  (.14)).  What  we  are  doing  here  is 
looking  at  the  ensemble  of  all  possible  realizations  of  N(t),  such  as 
Figure  2,  and  at  every  t selecting  out  all  those  which  have  had  an 
event  occur  in  (t , t + dt).  We  then  partition  all  theae  realisations  into 
three  subsets : 

(1)  those  with  no  hits, 

(2 ) those  with  only  one  hit,  and 

(3)  those  which  have  terminated. 

We  are,  of  course,  looking  at  subsets  of  (l),  (2),  and  (3)  where  N ■ n 

for  all  possible  values  of  N.  With  this  understanding  of  the  nature  of 

h°(t),  h*(t),  and  h (t)  we  may  now  write  three  interconnected  integral 

n n n 

equations.  For  example, 


h^(t )dt 

D 


dt  q [ f (x ) h°  (t  - x)dx 

J r\  n-1 
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which  states  that  if  do  hits  had  been  scored  in  n-1  rounds,  the  last  of 
which  was  fired  at  time  t - x,  then  if  we  fire  another  round  at  t it 
muot  fail,  in  order  ter.'  the  system  to  be  in  the  state  of  having  just  fired 
a round  which  was  the  n**1  one,  and  no  hits  have  been  made  yet.  The  inte- 
gral is  just  the  convolution  of  h^_1(t)  with  f(t).  If  we  now  define 

h°(t)  ■ ^ h£(t) 

n»l 

m 

■ q Jo  tix)  I Vx(t  ‘ x>ax 

n*l 

- qf(t)  + 4 ( f(x)  h°(t  - x)dx  . (21 ) 

J0 

Mote  that  for  n ■ 1 we  merely  have  one  selection  from  f(t),  multiplied 
by  the  probability  of  a failure,  in  order  to  get  h£(t)  (see  Equation  (l8)). 

In  a similar  manner, 


h*(t )dt 
n 

- dt  p J f(x)  h°_1(t-x)dx+  dt  q J f(x)  h^-1(t-x)dx  , 

where  we  have  accounted  for  two  mutually  exclusive  ways  to  get  to  the  state 
of  having  fired  one  round.  These  are,  no  successful  rounds  in  n - 1 
firings  and  a hit  on  the  nth,  and  one  hit  somewhere  in  the  first  n - 1 
rounds  and  no  hit  on  the  n^  round.  Sunning  as  before, 
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(22) 


fa1  (t ) 

■ pf(t)  + p [ f(x)  h°(t  - x)dx  ♦ q I fUJh^Ofc  • x)dx  , 
J0  J0 


where  we  note  that 


h*(t)  - 0; 
o 


and  finally, 


h (t  )dt 
n 


f(x)  h*  . (t  - x)dx 

U"1 


9 


which,  when  sunned  gives, 


h(t ) 


f(x)  hX(t  - x)dx 


(25) 


Equations  (2l),  (22)  and  (25)  nay  he  solved  in  that  order,  as  they  depend  on 
each  other  in  that  order.  Taking  the  characteristic  function  of  both  sides 
of  (21 ) and  denoting  the  cf  of  h°(t)  by  *°(u), 

♦°(u)  - q«P(u)  + q?(u)  *°(u)  , 


which,  when  solved  for  *°(u)  is 


,°(u)  . isSal 


l - qcp(u) 


(2k) 


Mow,  taking  the  cf  of  (22)  and  describing  the  cf  of  h^(t)  as  ^(u), 


^(u)  ■ pf(u)  + pf(u)  *(u)  + q«p(u)  ^(u) 
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substituting  #°(u)  fra n (2U)  and  solving  for  ^(u),  on*  obtains: 

*(u)  - • 0 

[1  - 4 ?(u)] 

Finally « taking  the  cf  of  both  sides  of  (23) 

*(u)  ■ p«P(u)  ^(u) 

and  substituting  ^(u)  from  (25)  into  this, 

•<«>  • [ 1^7  J ' ( 

which  again  is  the  same  as  (l6),  as  it  shcxild  be. 

The  advanta ge  of  this  technique  is  that,  in  spite  of  the  lengthy 
exposition  given  above,  it  is  often  the  simplest  and  fastest  way  to  a 
given  desired  result.  Also,  note  that  we  have  only  used  Equations  (l8), 
(19)  and  (20 ) for  definitions.  We  could  have  proceeded  directly  to  the 
solution  with  (20 ),  as  it  is  Just  the  mixture  solution.  Also,  we  can  use 
Equations  (2k)  and  (25)  (or,  (l8)  and  (19))  to  get  additional  information 
on  the  process.  These  have  physical  meanings  such  as:  the  inversion  of 
(2U)  will  give  the  (improper)  density  function  of  the  time  to  aero  hits, 
and  (2  5)  gives  the  pdf  of  the  time  to  one  hit.  These  additional  rela- 
tions may  also  be  obtained  by  the  mixture  technique,  but  each  would 
require  a separate,  Independent  calculation,  in  e»ctly  the  same  manner 
as  described  for  h(t)  in  Section  A above. 


D.  Renewal  Theory  Integro- Differential -Difference  Equation 
Techniques 

For  the  renewal  process  described  in  Section  C shore,  the  integro- 
dlfferential-difference  technique  nay  be  applied  if  the  interfiring  tines 
have  a negative  exponential  pdf  (ned).  The  reason  for  this  is  that  the 
ned  FV  has  no  memory  and  the  process  is  reduced  to  a terminating  semi* 
Markov  renewal  process  as  a result. 

However,  if  the  supplementary  variable  technique  (see,  Cox,  1955  or 
Keilson  and  Kooharian,  I960)  is  applied,  the  restriction  to  ned  IFT's 
is  removed  and  general  IFT's  may  be  considered  at  the  cost  of  some 
complication  in  the  mathematics.  In  this  fas.  Ion,  the  process  goes  from 
non- Markov  to  semi-Markov.  This  procedure  is  widely  used  in  the  literature 
pertaining  to  Stochastic  Duels,  but  much  ambiguity,  many  notatlonal  diffi- 
culties and,  in  came  cases,  unnecessary  complications  are  widespread. 
Consequently,  we  shall  illustrate  the  method  for  our  example  in  some  detail, 
in  the  hope  that  the  reader  will  find  the  literature  more  readily 
accessible. 

As  before,  the  process  may  be  considered  graphically  as  shown  in 
Figure  5.  The  situation  depicted  there  is  for  the  case  where  the  process 
is  observed  at  some  arbitrary  time,  t,  and  its  state  is  noted.  In  this 
case,  the  state  is  that  V(t)  - n and  there  have  been  no  hits.  The  time 
since  the  last  firing  is  a RV  and  is  in  the  interval  (t  - y , t - y + dy). 
The  corresponding  new  random  variable  Tn,  is  the  supplementary  variable 
and  in  renewal  theory  terms  is  the  backward  recurrence  time  (see,  Cox  and 
Miller,  1965,  p*  559).  The  other  two  cases  are: 

(l)  the  seme  as  shown  in  Figure  5,  except  that  one  hit  has  oc cured 


Backward  Recurrence 


At. 

on  one  of  the  n firings  (possibly,  even  on  the  n firing), 
and 

(2 ) at  the  time  of  observation,  the  process  has  terminated,  i.e., 
the  second  hit  occurred  on  the  n**1  round. 

In  this  latter  case,  the  supplementary  variable  i_  unnecessary  since  no 
further  events  can  take  place.  The  supplementary  variable  is  subscripted 
with  n as  it  is  a function  of  n,  but  its  realization  is  not  sub- 
scripted, since  ve  wish  to  look  backwards  the  same  distance  y for  all 
In's. 

The  properties  of  the  supplementary  variable,  Tn,  and  its  rela- 
tionships to  the  DT,  X (which  is  most  important),  are  veil  known  and  may 
be  summarized  as  follows : 


(a)  Tn  is  measured  to  the  left  from  t and  f„  (y  ;t)dy  is  the 

th  n 

probability  that  the  n firing  lies  in  the  interval 
(t  - y , t -y+  dy); 


(b)  fy  (y  ; t) 
n 


h^*(t  - y)  Fj(y),  0 < y < t. 


1 » 


is  the  already  familiar  pdf  of  the  renewal  process  for 
either  zero  hits  or  one  prior  hit,  with  N(t)  ■ n; 

(c)  P[y  < XBfl  < y ♦ dy  | Xnfl  > y]  - P[y  < X < y ♦ dy  | X > y] 


fx(y)«Jy 

- x(y)dy  - ■ ■ , 

p£(y) 


which  implies  that  the  probability  of  a firing  in  the  interval 

(t , t + A)  - X(y)A  + o(A),  where  lim  ■ 0 ; 

A -+0  J 


i and 


...  . 


I 


z 


% 


-£  x(t)di 

(d)  fx(y)  ■ X(y)  e 

(e)  from  (c)  and  (d)  above,  F°(y) 


-Jjj  X(Od6 


With  this  background,  ve  now  define  a function  which  is  the  probability 
that,  if  the  process  is  observed  at  some  arbitrary  time  t,  with  n 
firings  having  occurred,  and  no  hits  have  been  made,  and  that  the  last 
firing  was  between  y and  y-dy  time  units  earlier. 


H°(t , y)dy  - P[N°(t)  «n  , y < Y&  < y - dy] 

j P[T®  < t , y < Tn  < y - dy] ; 0<y<t,  n>0 

I 0 ; y > t > 0,  n>0 

H°(t  , y)dy  - 6(y  - t)  F^(y)  ; y - t > 0 . (2T) 

The  first  form  of  (27)  is  to  be  interpreted  as  a Joint  probability  mass 
function  (improper)  on  the  number  of  romds  fired  with  no  hits  occurring 
(N°(t)),  and  a probability  density  function  on  the  time  since  the  last 
firing,  Yn,  with  a parameter  t.  The  second  (equivalent)  form  is  inter- 
preted as  a Jo'nt  distribution  function  (improper)  on  T°,  the  time  of 

n 

firing  the  n^*1  round  with  no  hits  and  a pdf  on  Yn,  with  a parameter 
n.  We  also  note  that  when  n ■ 0,  that  y must  equal  t,  which 
accounts  for  the  Dirac  delta  function  6(y  - t). 

In  an  exactly  similar  fashion  we  have, 
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(2 


H*(t  , y)dy  - PtN^t)  «n  , y < < y-  dy] 

P[lJ  <t,y<Yn<y-dy];  0 < y < t,  n > 1 

9 

0 » y > t,  n > 1 

where  N^t)  Is  the  RV  number  of  rounds  fired  with  one  hit  somewhere, 
and  where  is  the  RV  time  of  firing  of  the  nth  round  with  one  hit 

in  one  of  the  n rounds  fired.  Now,  the  case  where  the  process  has 
terminated  is  described  by 

H (t)  * P[N(t)*n,  2nd  hit  on  the  nth  round] 

n 

■ P[Tn  < t , process  terminated] ; t > 0,  n > 2 . (2 

In  this  case,  the  first  form  is  a proper  pof  on  N with  parameter  t and 
the  second  form  is  a proper  df  on  In  where  n is  a parameter. 

The  fact  that  certain  of  the  mass  and  density  functions  above  are 
improper  will  be  seen  as  we  derive  them.  As  in  Section  C,  above,  this 
occurs  whenever  we  consider  a subset  of  the  sample  space  which  is  not  t1ust 
the  terminated  subset.  The  reason  for  this  is  not  easy  to  see,  but  it  is 
related  to  the  fact  that  all  realizations  terminate  with  probability  one. 

With  this  background  we  may  now  derive  the  integro- differential- 
difference  equations  which  govern  the  process.  If  we  ask  about  the  state 
at  some  forward  time,  t+£,  note  that  this  also  simultaneously  extends 
y to  y+  A. 

Now,  from  elementary  calculus  and  referring  to  Figure  4,  if  we  ask 


| 

i » 


-2 


» 


Hn(t+A,  y+A) 


Hn(t,  y)  Hn(t,  y+A) 


FIGURE  k 

about  being  in  stete : n firings  and  no  hits  at  t+A,  we  have 

H°(t  + A , y + A) 

■ H^(t  , y)  + -^  H^(t  , y)A  + H°(t  , y+  A)A  . 

Also,  this  same  probability  is  that  we  are  in  the  state:  n firings  and 
no  hits  at  t,  and  no  firing  occurs  in  the  interval  (t , t+A) , i.e., 

H°(t , y ) [1  - X(y)A].  Equating  these  two  statements,  we  have, 

H°(t  , y)[l  - X(y)A] 

* H°(t , y)  + H°(t  , y)A  + H°(t , y+  A)A  . 
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y Rearranging  terms , dividing  through  by  A,  and  taking  the  limit  as  A -»  0 

gives, 

^ + -Jf  + X(y)  J H°(t,y)«0;  n-  0,1,2,...  , 0 < y < t.  (30 ) 

Kote  that,  although  H°(t,  , y)  is  quite  different  in  form  from  H°(t , y) 
for  n > 0,  it  still  satisfies  (50 ). 

In  a precisely  similar  manner,  ve  can  ismediately  write  down  for  the 
state:  n firings  and.  one  hit, 

( + + X(y)  J H*(t  ,y)  - 0;  n - 1,2, . . . , 0 < y < t . (3l) 

For  the  situation  of  n firings  and  two  hits  (process  terminated),  there 
are  two  mutually  exclusive  possibilities  when  we  go  to  t+  A from  t. 

Either  the  process  had  terminated  before  t (and  therefore,  remains 
terminated),  or  it  had  not  terminated  earlier  and  does  terminate  in 
(t,t+A),  thus, 

H (t  + A)  » H (t)  + [ hJ  ,(t  , y)  X(y)Apdy  ; 

n n j q n-i. 

n > 2,  t > 0 . 

■«  - “ “ 

i 

{ Rearranging  terms,  dividing  through  by  A and  taking  the  limit  as  A -♦  0, 

•# 

-1 

' 5;  t 

[Hn(t)l  * hn(t)  -p  J ^.i(t  i y)  X(y)dy  ; 

n > 2,  t > 0 . (32  ) 

S 
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since  there  is  no  way  to  have  a hit  without  a firing.  The  two  terns  on  the 
right-hand  side  of  (37)  account  for  the  two  ways  to  get  to  one  hit  at 
exactly  tine  t,  i.e.,  no  prior  hits  and  get  a hit,  and  secondly,  have  a 
prior  hit  and  get  a miss  at  t. 

We  now  define  the  following  functions 


H°(t , y) 

- 2 
n«0 

<(t,y) 

H^t  , y) 

«• 

- 2 

«J(t , y) 

n*l 

h(t) 

m 

- 2 

h (t ) 

T& 

n 

h°(t) 

« 

• 2 

h°(t) 

n»l 

n 

hX(t ) 

■ 

s 

^(t)  j 

Mote  that  upon  stunning  on  n, 

H°(t , y )dy  « P[T°  <t,y<T<y-dy]  , 

where  T°  is  the  RV,  time  to  a firing  (for  any  n)  with  no  hits  and  T 
is  the  RV,  time  since  last  firing  (for  any  n).  Similarly,  H*(t,y) 
involves  T1  the  RV,  time  to  a firing  (for  any  n)  with  one  prior  hit. 


Another  important  point  is  that 


11m  H°(t , y)dt  - P[T°  <t,0<T<0-dy] 
y -»0 

- P[T  < T°  < t + dt] 

- 6(-t)dt  + h°(t )dt  , 


from  Equations  (27)  and  (39  )•  This  Is  because  the  definition  of  11° (t , y) 
Includes  the  case  of  n ■ 0 and  h°(t ) does  not.  However, 


11m  H^tjy)  - h*(t)  . 

y -»0 


Vow,  upon  summing  on  n and  applying  definitions  (39)  to  Equations  (30 ), 
(31),  (32),  (35)  and  (37), 


( + ■Jji  + ) h0^  * * °*  0 < y < t, 


( ^ + 'k  + * 0 # 0 < y < t, 


h(t)  - p / H1(t,y)  X(y)dy  , t > 0, 

J0 


h°(t ) - q f H°(t , y)  X(y)dy  , t > 0, 

J0 


hX(t)  7 q [ H^t , y)  X(y)dy  + p [ H°(t , y)  X(y)dy  , t > 0 , 
J 0 J0 


end  finally,  the  Initial  condition  from  (27), 


This  cooes  about  because,  when  t ■ 0,  y ■ 0 and  F^(o)  • 1. 

Before  proceeding,  it  must  be  noted  that,  in  aost  of  the  literature, 
the  preceding  step  is  not  a simple  svumnation  on  n,  but  '■ether,  the  form- 
ation of  a geometric  transformation  (sometimes  referred  to  as  a probability 
generating  function  or  Z transform)  first*  This  is  done,  by  say, 

z11  H^(t , y)  where  a is  the  transform  variable.  In  almost  all  cases, 
this  step  is  unnecessary  and  yields  nothing,  since  authors  usually  set 
Z ■ 1,  and  we  are  right  where  we  are  now.  The  only  reason  for  using  the  Z 
transform  technique  at  this  point  is  if  one  wishes  to  discover  some  of  the 
properties  of  the  peif  on  N(t).  This  may  be  done  in  the  usual  fashion  by 
taking  appropriate  derivatives  of  the  Z transform  with  respect  to  s . 


The  next  step  is  to  convert  all  functions  to  their  characteristic 
functions,  in  the  variable  t,  by  the  appropriate  transformation.  We 
define : 


cf  *ls : 


H°(t , y) 

~ 7°(u,y) 

HX(t , y) 

~ TX(u , y) 

h(t ) 

~ ®(u) 

h°(t) 

~ *°(u) 

hX(t ) 

~ ^(u) 

0 


In  performing  thc*?e  transforms,  we  note  that  the 

cf  ( a*H°(t,y)J  - - iu  ^(u  , y)  - H°(0  , y)  . 


Proceeding  to  operate  on  Equations  (42 ) through  (46)  and  using  (47), 


y > o , 


f - iu  X(y)  ) *°(u  , y)  » a(y)  , 

^ - lu  ♦ X(y ) J ^(u  , i ) » 0 , 


y >o  , 


*(u) 

♦°(u) 


^(u) 


" P [ ?1(u,  y)  X(y)dy 

J0 

■ q [ *°(u,  y)  X(y)dy 

J0 

q [ ^(u  , y)  X(y)dy 
J0 

> f *°(u  , y)  X(y)dy 

A 


♦ P 


I 


> 


We  shall  Illustrate  how  the  right-hand  side  (rhs)  of  the  last  three  equa- 
tions come  about  by  looking  at  the  derivation  of  the  rhs  of  (5l). 


rhs 


iut 


dt  f H^(t , y)  X(y)dy 

w n 


dt  J H^Ct , y)eiut  X(y)dy 


low,  since  y < t,  the  integration  above  is  over  the  shaded  region 
depicted  in  Figure  5>  first  in  the  y direction  and  then  in  the  t 
direction.  Reversing  the  order  of  integration  ve  have 


rhs 


P 


/ X(y)dy  [ ciut  H^t  , y)dt 


y 


FIGURE  5 

Integrating  to  Obtain  Equation  (5l) 

and  since  H^(t , y)  - 0 for  t < y,  the  inner  integral  is  ^(u  , y ) and 
ve  have  the  desired  result. 

The  problem  now  is  to  solve  the  simultaneous  set  of  equations  (49) 
through  (53).  Equation  (49)  is  a linear,  nonhomogeneous  differential  equa- 
tion in  the  variable  y.  Biis  may  be  solved  in  the  usual  way  by  getting  the 
general  solution  to  the  homogeneous  equation  (rhs  • 0)  and  adding  a par- 
ticular solution  for  the  entire  equation.  Thus,  setting  rhs  ■ 0, 


*£(»,y) 
d In  *°(u  , y) 

— fc — 


iu  - X(y)  , or 


iu  - X(y)  , or 


(54a) 


(54b) 


-32- 


In  »°(u,y) 

- la,  - J 

nr 

X(Odt  ♦ c , 

0 

(5^c) 

*?»•*>  ■ 

x - !l 

X(i)dS 

(5^) 

where  K end  e are  tone  constant!,  and  Y°  ic  the  general  aolvtlon  to 
the  homogeneous  equation. 

low,  let  us  try  the  following  particular  solution,  Y° , 

T°(u,y>  - iKy).1”*  - 'l  X(S)dS  (55) 

where  U(y)  is  the  unit  step  function,  in  (^9),  and  we  see  at  once  that  it 
is  satisfied.  Therefore, 

*°(u,y)  - - (U(y)+  K]eiuy“^o  X(t)dfe  . (56) 

We  notice,  from  (Uo)  and  (1*6),  that  if  we  let  y -*  0, 

Y°(u)  - 1 ♦ *°(u)  , 


which  immediately  gives  us  that  K ■ *°(u).  Thus, 


*°(u,y) 


[U(y)  + t°(u)]e 


iuy  - f*  X(6)di 


(57) 


The  solution  to  (50)  follows  iamediately,  since  there  is  no  complicating 
rhs  and  we  nay  write 


(58) 


•-t1(u  , y) 


Now,  aubstituting  (57)  into  (52), 

r-  _ iuy - / y X(i)d£ 

♦°(u)  - q / [u(y ) + •‘’(uMe  0 X(y)dy  . (59) 

J0 

Remembering  that  U(y)  - 1 for  y > 0,  the  integral  an  the  rhs  is  Just 
the  cf  of  the  pdf  of  X.  Therefore,  (59)  is 

*°(u)  - q(l  + *°(u)3  «p(u)  , 


from  which 


•°(u) 


life) 

l - q?(uj 


(60) 


Substituting  (60)  into  (57) 


70(u,y) 


U(y) 


a ?(u) 


iuy  - /0y  X(fc)d£ 


(61) 


Substituting  Iquations  (58)  and  (6l)  into  (55) 


X(y)dy 


^(u)  - q**(u)  J 


- iuy-/yX(Odfc 
e ° 


♦ P 


I"  [«<*>♦  r*" 


iuy-  x(Odt 


or 


^(u) 


q«^(u)  qp(u)  + pf(u)  1 + 


i -q«p(u; 


or 


^(u) 


Pffr) , 

II  -<lf(u)]2 


which,  for  (58)  gives 


ll-qvMf 


Finally,  (63)  is  substituted  into  (5l) 


♦(„)  - ! [ . 

[1  - q?(u)]  J0 


m iuy-jfyX(6)dt 


X(y)dy 


or 


•(u) 


.,».*.(»  L 

1-qfIuT 


t 


which  1b  the  desired  result  and  checks  with  all  previous  derivations. 

At  this  point  i one  might  question  the  need  to  use  such  an  involved 
and  complicated  method  to  get  the  result  in  Equation  (6*0  when  other 
methods  are  simpler.  The  answer  to  this  is  that  if  it  is  the  only  result 
desired,  then  there  is  no  advantage.  However,  it  should  be  noted  that 
the  reward  for  the  additional  effort  is,  as  it  should  be,  that  additional 
information  is  available,  namely,  ore  can  invert  Equations  (60),  (6l),  (62) 
and  (63)  to  get  probability  functions  on  the  distribution  of  times  for 
zero  hits,  one  hit,  joint  times  for  zero  hits  and  time  since  last  hit,  and 
one  hit  and  time  since  last  hit.  Also,  as  previously  noted,  if  the 
geometric  transform  is  used,  information  on  the  number  of  rounds  fired  can 
be  obtained. 
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APPENDIX 


SOME  USEFUL  RESULTS  IN  THE  THEORY  OF  CHARACTERISTIC 

FUNCTIONS 


0.  Introduction 

This  appendix  contains  sane  theorems  and  other  results  from  the 
theory  of  characteristic  functions  which  are  useful  in  derivations  or 
applications  in  the  Theory  rf  Stochastic  Duels.  All  of  these  results, 
except  Theorem  (A15),  are  from  a compilation  in  bibliographic  item  A6. 

The  proof  of  Theorem  (A15)  is  given  herein. 

We  use  ^(u)  to  denote  the  characteristic  function  of  the  pdf 
fg(t).  Except  for  the  lfcrseval  theorems,  the  results  are  for  positive 
RV's  only.  The  notations  /L  and  are  the  same  as  those  given 
earlier  in  this  work. 

1.  Some  Bsrseval  Theorems 

Three  useful  versions  of  Ihrseval's  Theorem  are  given  below. 

(Al) 

(A2) 

(A3) 


/.  fA(t)  fB(t)dt  * £i  j[.  Vtt)  Vu)du 

/.  •tak  rkM  fB(t)dt  mh  I Vu-’,)  hMi"  ■ 
J elut  fA(t)  fB(t)  fc(t)dt 

' *7  l.  ,*<u'v)(  l.  Vv*t)  ’’o(t)4v)*'- 


A-l 


2.  Properties  of  Characteristic  Functions  of  Positive 
Fandom  Variables 


In  what  follows,  only  positive  random  variables  are  considered, 
i.e.,  pdf's  such  that 


and 


f(t)  >0,  t > 0 \ 

f(t)  - 0,  t < 0 J 


f (t  )dt 


1 


with  characteristic  function  and  inverse 


?(u) 

f(t) 


" J *iUt 

& C • ** 


real. 


(AUa) 


(AUb) 


(A5) 


respectively. 

The  properties  of  interest  follow. 

?(0)  - 1 , (A 6) 

|f(u)|  < 1,  Imaginary  u > 0 . (A7) 


This  Implies  no  singularities  in  the  upper  half  of  the  complex  plane. 


l?(-u)|  < It 


(A 8) 


This  Implies  no  slngulsritles  In  the  lower  helf  of  the  complex  plane. 


|?(u)|  < -r-r  , f(t)  a differentiable  function  (A9) 

1 1 of  bounded  variation.  Imaginary 
u > 0,  k ■ posit ire  constant. 


This  implies  that  f(u)  diminishes  as  l/R  in  the  upper  half-plane  where 
R is  the  radius  of  a semicircular  path  of  integration  in  the  complex 
plane. 


k(-u)|  < -I — r , f(t)  a differentiable  function  (Aio) 

*u'  of  bounded  variation.  Imaginary 
u < 0,  k ■ positive  constant. 


This  implies  that  f(-u)  diminishes  as  l/R  in  the  lower  half-plane,  R 
the  same  as  in  (A9). 


! Imaginary  u > 0 

(AU) 

Teaginary  v < 0 


This  implies  no  singularities  in  the  upper  half  of  the  complex  u plane 
and  the  lower  half  of  the  complex  w plane. 


|*(u  - v)  1 < -i r , f (t ) a differentiable  function  (A12  ) 

•u  " of  bounded  variation.  Imaginary 
u > 0,  imaginary  w < 0,  and 
k ■ positive  eanstantT 


This  implies  that  f (u  - v)  diminishes  as  l/R  in  the  upper  half  of  the 
u plane  and  in  the  lower  half  of  the  v plane.  R,  as  shove. 
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3.  Bane  Theorems  Involving  Character  la  tic  functions  of  fbsitlve 


id  a Variables 


We  hare  that 


r *iut  **  - a jTi  ■ ■ ■ 


-iva , 


dv 


(A13) 


Rote  that  for  u ■ 0,  this  also  gives  an  express  ion  far  the  distribution 
function  of  a random  variable  In  terme  of  characteristic  functions.  Also, 


r .ittt f(t)dt ■ *<u) . ^ r - «~i”) *, . (Ait) 

^ a - “ 


Rote  that  far  u - 0,  this  provides  an  expression  for  the  compleaentary 
distribution  function  In  terns  of  characteristic  functions. 

The  characteristic  function  of  the  distribution  function  is  given 


( J*  f U)d  t j dt  - - ^ , Imaginary  u > 0 


(A15) 


ftroof : integrating  by  parts, 


a: 


fU)dt  )dt 


iut  rt  i* 

■v/0  ^H.0- 


r " iut 

Jo  V f(t)4t 


y(u) 

• i..  • 


low,  to  explore  the  first  tern,  let  u be  analytically  continued  into  tbe 
complex  plane.  Sms,  u ■ x + ly,  x,y  real. 


i 

iu 


11a 

t-** 


iut 


■jL  lia  e1**  e"3*  . 


t -» * 


Since  le1^!  ■ 1, 


15  lim  #iUt 


0 , y > 0,  any  x 

ixdetennlnate,  y ■ 0,  any  x 
* , y < 0,  any  x 


Continuing  our  listing  of  results. 


Q.E.D. 


/o’  etat(/t’  tU)il)at  • ^ 


(A16) 


This  is  the  characteristic  function  of  the  complementary  distribution 
function. 


i: 


fB<l)4!)dt  - J -a 


• f.(u-v)[*L(v)-l] 

— r2 d» , (ait) 


/0"v*>  * - A /]  - 


A(-u)[t>  (u)  - 1] 

A s du  , (A18) 
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A S 4u  , (A2Z) 

w 


vCiTTi.' 


j(l)4l  J dt 


A.  r 

m Jm  m 


-iub  - / V [*B(u)  ’ 11 

• • - t; 


du  i 


fB(t)dt  Jdt 


(£ 


r _ y ✓ f-  fA(u-aw)  •“ttw  [^(w)-l] 


dv  Jdu 


j(e)di  J 


fA(-»u)[<PB(u)  - 1] 


ut  - ii  f Fm  Ta^u’  - 1] 


*u: 


dw  Jdu 


•(,)dl  X /»"  fC(,,)  [ X fD(p)ie]  *•  j 


• t«(-u- ▼){?«(▼)  - lldv 


?.(u  - *)[?„(*)  - lldw 


Cur  final  result  la  that  for  any  Integral  of  the  fora 


P[A] 


tA(-u)  t^(u) 
u 


du 


> 


where  there  nay  ha  may  nunber  of  of 'a  of  positive  RV'a  la  the  integrand, 
as  long  as  at  least  oaa  baa  a poaltiva  argument  and  at  laaat  one  ha;;  a 
negative  argument  and  all  are  divided  by  u,  then 

P[A]  ■ 2 f Hi  Integrand  (A27) 


where  C la  any  path  la  the  lower  right  half  of  the  complex  u plane, 
which  atarta  at  u - 0 and  terminates  at  + * • The  path  auat  remain  on 
the  real  u aria  or  be  between  the  real  u axle  and  the  nearest  singular- 
ity la  the  lower  half  plane.  Examples  are  shown  la  Figure  Al. 


FIGURE  Al 
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PART  II 


ANNOTATED  BIBLIOGRAPHY  OF 
RESEARCH  PAPERS 


This  part  of  the  work  contains  a comprehensive,  fully  annotated 
bibliography  of  all  research  results  known  to  the  author.  An  exhaust- 
ive search  of  the  free-world,  literature  has  been  made.  Some  manuscripts 
which  have  never  been  published  in  the  open  literature  are  also 
included . 

A few  papers  included  deal  primarily  with  many  - versus  - many 
situations,  but  contain  some  material  on  one  - versus  - one  duels.  In 
this  case,  only  the  latter  is  annotated  and  appears  herein. 

The  format  is  as  follows : 

(l)  Each  paper  has  an  alpha-numeric  designation  consisting  of 
the  first  one  (two  if  necessary)  letters  of  all  authors 
last  names,  plus  a number. 

The  papers  are  arranged  alphabetically  by  first  author,  starting  with 
Bole  authorship  and  then  dual  authorship,  and  so  on.  The  number  indi- 
cates the  chronological  sequence  of  all  papers  by  a specified  author, 
or  set  of  authors.  The  chronology  is  based  on  the  date  of  the  first 
citation  given  for  each  entry  in  the  bibliography.  There  may  be  more 
than  one  citation  given,  as  it  is  customary  to  produce  first  a company 
or  university,  etc.,  version  and  then  an  open  literature  publication. 
The  latter  is  always  the  first  citation.  The  reason  all  versions  are 
given  is  that  a reader  who  may  have  access  to  one  version,  may  rest 
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assured  that  there  Is  no  difference  in  content  between  versions,  except 
as  noted. 

Next,  if  there  are  both  fundamental  marksman  (Fm)  and 
fundamental  duel  (FD)  results,  the  FM  part  is  given  first  and  the 
FD  portion  is  next.  The  format  for  each  part  is  the  same,  and  goes 
as  follows: 

(1)  Identified  as  FM  or  FD 

(2 ) For  each  model  considered  under  FM  or  FD 

(a)  identified  as  CRIFT  or  FIFT,  or  mixed 

(b)  all  modifications,  e.g.,  ammunition  limitations,  etc. 

(3)  General  solutions  which  are  given 

(4)  ftirticular  examples  with  details 

Finally,  at  the  bottom  left-hand  side,  a listing  of  all  the  prin- 
cipal techniques  of  derivation  used  is  given. 


The  fundamental  marksman  (FM)  problem  is  defined  as  a marksman 
firing  at  a passive  target.  He  fires  at  certain  intervals  and  either  hits 
or  misses  on  each  round.  The  trials  cease  on  the  first  hit.  The  hit 
probability  is  constant  from  round  to  round.  The  time  between  rounds  may 
be  a continuous  RV  or  a constant.  He  starts  with  an  unloaded  weapon, 
has  unlimited  time  to  hit  the  target  and  has  unlimited  ammunition. 

The  fundamental  duel  (FD)  pits  2 marksmen,  as  defined  above, 
against  each  other.  The  duel  terminates  when  one  hits  the  other,  or  both 
are  hit  simultaneously.  They  start  at  the  same  time  and,  in  general,  have 
different  kill  probabilities  and  different  interfiring  times. 

NOTATION  AND  ABBREVIATIONS 


CRIFT  - continuous  random  interfiring  times 

cf  - characteristic  function  ■ " eilVt  f(t)dt 

CRV  - continuous  random  variable 

df  - distribution  function 

DDC  - defense  documentation  center  (documents  may  be  ordered  from 

this  center  using  numbers  as  shewn  in  bibliography) 

Erlang  (k)  - a RV  with  pdf  given  by 

v k-1 

fX(x)  m 'Gr-“l)'{  e ax*  x > 0,  k - 1,2,... 

■ 0 , elsewhere 

E(X]  - expectation  of  X 


B-3 


• fundamental  duel 


FD 

F1FT  - fixed  interfiring  times 

FM  - fundamental  marksman 

gt  - geometric  transform  - 2 f (n)z11 , also  sometimes  called  the 

z transform  and  if  f(n)  are  elements  of  a pmf,  sometimes 
called  a probability  generating  function 

h(t)  - pdf  of  RV,  marksman's  time  to  a kill 

H(t)  - df  of  RV,  marksman's  time  to  a kill 

H - the  event  of  a hit 

H - the  event  of  no  hit 

IFT  - interfirlng  times 

ild  - independent,  identically  distributed 

IT  - Laplace  transform  ■ e-st  f(t)dt 

mgf  - moment  generating  function  (same  as  the  Laplace  transform 

with  s replaced  by  -s) 

MLE  - maximum  likelihood  estimate 

ned  - negative  exponential  pdf,  given  by 

fx(x)  - r e”rx  , x > 0 

■ 0 , elsewhere 

K - number  of  rounds  fired  j may  or  nay  not  be  a RV  (sometimes 

used  as  a general  constant,  as  noted  in  text) 

2 o 

V ( u , c ) - normally  distributed  RV  with  mean  u and  variance  cr" 

pdf  - probability  density  function,  denoted  by  fx(x) 

pmf  - probability  mass  function,  denoted  by  p^x) 

pgf  - probability  generating  function  (or,  sometimes  called 

geometric  transform  or  z transform) 

P[A]  - probability  A wins  the  duel 


B-ifr 


I 


P[B]  - probability  B wins  the  duel 

P[AB]  - probability  neither  A nor  B wins  the  duel,  consequently, 
a draw 

r - rate  of  fire  or  value  of  RV,  R 

P,  - RV,  number  of  hits  to  a kill  (for  a situation  where  more 

than  one  hit  is  required) 

PV  - random  variable 

T - PV,  marksman's  time  to  a hit  (kill)  - may  be  subscripted  as 

appropriate 

TCSF  - time -of- flight  (may  be  a RV) 

V[X]  - variance  of  X 

X - interfiring  time  RV 


! 

i • * 

V' 

I ; * 

S^jfrrp^:: 
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Al.  AMCKER,  C.J.,  Jr.,  "Stochastic  Duels  With  Limited  Aromunition  Supply," 
Operations  Research,  Vol.  12,  Mo.  1,  Jen. -Feb.,  1964,  pp.  36-50. 

(Also,  System  Development  Corporation,  Sente  Monica,  CA.,  Document 
SP-1017/001/00,  23  April,  1963,  24  pp.,  DDC  Mo.  AD-404  104.) 

FD  - CRIFT  - (l)  Random  ammunition  supply 

- (2)  Fixed  ammunition  supply 

- (3)  Special  ease;  seme  as  (l)  above,  except  duel  terminates 

when  either  side  runs  out  of  ammunition  before  a kill 

General  Solutions:  P[A],  P[B),  P[AB]  (both  sides  run  out  of 

ammunition,  (l)  and  (2);  either  side  runs  out  of 
ammunition,  (3)  ) 

Examples : (l)  and  (2)  only 


Distribution  of  Humber  of  Rounds 


A 

Geometric 

Poisson 

Binomial 

Geometric 

Fixed 


B 

Geometric 

Geometric 

Geometric 

Geometric 

Fixed 


IFT  pdf's 

ned 

ned 

ned 

Erlang  (2) 
ned 


Curves  comparing  this  duel  with  FD 


mixture  technique 
characteristic  functions 


A2.  ANCKER,  C.J.,  Jr.,  "Stochastic  Duels  of  Limited  Time-Duration,"  CQRS 
Journal  (Canada),  Vol.  4,  Mo.  2,  July,  1966,  pp.  69-81. 

(Also,  System  Development  Corporation,  Santa  Monica,  CA.,  Document 
SP-101T/005/00,  30  March,  1964,  24  pp.,  DDC  Mo.  AD-436  529.) 


FD  - (l)  CRIFT 
- (2  ) FIFT 


Length  of  Duel 

(a)  Continuous  random  time  limitation 

(b)  Fixed  time  limitation 

(a)  Continuous  random  time  limitation 

(b)  Fixed  time  limitation 
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A5. 


X 


t 


(cont’d) 

General  Solutions : P[A],  P[Bl,  P[AB]  (both  sides  run  out  of  time 

or  kill  simultaneously) 

Examples : 

Distribution  of 


A’s  IFT 

B's  IFT 

Time  Limitation 

ned 

ned 

ned 

Erlang  (2) 

Erlang  (2) 

ned 

ned 

ned 

fixed 

Erlang  (2) 

Erlang  (2) 

fixed 

fixed 

fixed 

ned 

fixed 

fixed 

fixed 

Curves  comparing  this  duel  with  FD 


mixture  technique 
characteristic  functions 
number  theory 

ANCKER,  C.J.  Jr..  "Stochastic  Duels  With  Time-of-Flight  Included," 
OPSEARCH  (India),  Vol.  3,  No.  2,  1966,  pp.  71-92. 

Errata  OPSEARCH  (India),  Vol.  3,  No.  3,  1966,  p.  155. 

(Also,  System  Development  Corporation,  Santa  Mcnica,  CA. , Document 
SP-1017/009/00,  19  May,  1966,  28  pp. ) 

3 Procedures  Considered: 

(1)  No  Delay  - firing  proceeds  as  rapidly  as  possible,  no  delay  to 

observe  effect 

(2 ) Delay  - each  round  is  allowed  to  hit  before  next  round  is 

prepared  and  fired 

(3)  Mixed  - one  side  delays,  the  other  has  no  delay 

FM  - CRIFT  - (l ) No  Delay  - random  TOF 

- fixed  TGF 

- (2 ) Delay  - random  TOF 

- fixed  TOF 

FM-FIFT  - No  Delay  - fixed  TOF 

General  Solutions : (pdf  time  to  fire  killing  round,  pdf  time  to 

kill) 


t 
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A3,  (cont'd) 

Examples : 


Procedure 

No  Belay 
So  Delay 
Delay 


ITT  pdf 

ned 

ned 

ned 


TQF  pdf 


ned 

constant  (fixed) 
ned 


F 

FD  - CRIFT-  (1) 

Mo  Delay 

- random  TOF 

- fixed  TOF 

f 

i 

[ 

- (2  ) 

Delay 

- random  TOF 

i 

! 

- (3) 

Mixed 

- A fixed  TOF  (delay),  B 

- Both  fixed,  A deity,  B 

l 

h 

1 

FD  - FIFT  - (l) 

No  Delay 

- fixed  TOF 

- (2  ) 

Delay 

- fixed  TOF 

General  Solution:  P[A],  P[B],  P[AB]  (both  killed) 

Examples : 


Special  Case: 

FD  - CRXFT  - no  delay  - TOF  varies  linearly 
General  Solution:  P[A],  P[B]  P[AB] 
Example : 

A's  IFT  pdf 
ned 


B's  IFT  pdf 
ned 


mixture  technique 
characteristic  functions 
number  theory 


Procedure 

A's 

IFT  pdf 

A's 

TOF  pdf 

B’s 

IFT  pdf 

No  Delay 

ned 

ned 

ned 

ned  j 

No  Delay 

ned 

fixed 

ned 

fixed  j 

No  Delay 

fixed 

(c  x B's  IFT) 

fixed 

fixed 

fixed  I 

5 

Delay 

ned 

ned 

ned 

ned  | 

A Delay 

ned 

fixed 

ned 

aero 

A Delay 

B no  Delay 

ned 

fixed 

ned 

fixed 

| 

A4.  ANCKET,  C.J.,  Jr. , "Stochastic  Duels  With  Round-Dependent  Hit 

Probabilities ,"  Maval  Research  Logistics  Quarterly,  Vol,  22,  Mo.  3, 
Sept  , 1975,  PP-  575-583. 

(Also,  University  of  Southern  California,  Los  Angeles,  CA.,  ISE 
Department  Technical  Report  TR  74-3,  2 August,  1974,  16  pp.) 

FM-CRIFT  - Round  dependent  hit  probabilities 

General  Solution:  pdf  time  to  a kill 

FD  • CRIFT  - Round  dependent  hit  probabilities 

General  Solution:  P[A]  P[B] 

Examples : 


A's  IFT  pdf 

A’s  Hit 
Probability 
on  j - th  Round 

B's  IFT  pdf 

Hit 

Probability 

ned 

a -fM.  1 V,  qA  , 

ned 

fixed 

q J 1 J (n-i) 

(M,  a fixed  Integer) 

Erlang  (2 ) 

qA 

qj  T 

ned 

fixed 

ned 

a -*A 

qJ  J 

ned 

fixed 

Curves  for  last  example  and  compering  1-st  vith  FD 


mixture  technique 
characteristic  functions 


A5.  AMCKER,  C.J.,  Jr.,  "Stochastic  Duels  With  Bursts,"  Maval  Research 
Logistics  Quarterly,  Vol.  23,  Mo.  4,  Dec.,  1976,  pp.  703-711. 


(Also,  University  of  Southern  California,  Los  Angeles,  CA.,  ISE 
Department  Technical  Report  TR-73-5,  Nov.,  1973,  14  pp. ) 


I 


A5.  (cont‘d) 

FD  - Nixed,  CRIFT  - FIFT 

General  Solutions : P[A],  P[B] 

Example : A - FIFT,  B - ned  IFT 
Curves  of  solution 

FM-  (l)  Bursts  of  fixed  sise  V (random  time  between  rounds  and 
random  times  between  bursts) 

- (2)  Same  as  (l)  with  fixed  times  between  rounds  in  a burst 

General  Solutions  : pdf  time  to  a kill 

FD  - A fires  burst  of  fixed  size  N.  Rounds  within  a burst  equally 
spaced.  Random  time  between  bursts.  B is  CRIFT,  no  bursts. 

General  Solutions : P[A],  P[B] 


Example : 


A B IFT 

ned  between  bursts  ned 

Solution  curves,  comparison  with  FD 

mixture  technique 
characteristic  functions 


! $ 


A6.  AWCKER,  C.J.,  Jr.,  "Theory  of  Stochastic  Duels  - Miscellaneous  Results," 
TRASAHA  TECHNICAL  MEMORANDUM  2-77,  March,  1978,  39  PP*  > U.S.  Army  TRADOC 
Systems  Analysis  Activity,  White  Sands  Missile  Range,  New  Mexico,  DDC  No,. 
AO- 52158. 

FM  - Erlang  (n)  CRIFT 

General  Solution:  pdf  and  df  time  to  a kill 

Examj&e:  n ■ 1,  n • 2,  solution  curves 

FD  - Tactical  equity  (each  side  fires  first  l/?  the  time  and  then  FD 
starts) 

General  Solution:  P[A] 


Example : Erlang  (2)  CR.TT's 


A6.  (ccnt’d) 


FD- CRIFT's  - tactical  equity  (except  2-nd  fixer  returns  fire 
immediately) 

General  Solution:  P[A] 

Example : ned  IFT's 

FD  - CRIFT  - initial  surprise  a CRV 

Example : IFT's  - ned,  surprise  pdf  - Laplace 
- solution  curves,  P[A] 

Sub -Example : 

Surprise:  (l)  pdf  - ned 

(2 ) pdf  - ned  for  negative  time  only 
FD  - Erlang  (n  , m)  CR IFT's 

General  Solution:  P[A],  P[B] 


6 problems  by  Thompson  in  H are  simplified  by  using  characteristic 
functions,  viz: 

(l)  (a)  (b)  (c)  and  (2)  (a),(b),(c) 

1 problem  by  Thompson  in  12  is  simplified  by  using  geometric  transforms 
and  characteristic  functions,  viz:  2 (c) 

Some  useful  results  in  the  Theory  of  Characteristic  Functions  are 
listed  or  derived  as  follows : 

3 Parseval  Theorems 

7 properties  of  cf's  of  positive  RV's 
15  theorems  concerning  cf's  of  positive  RV's 
2 theorems  on  contour  integration  in  the  complex  plane  for 
special  integrands  are  developed*  Useful  in  numerical 
integration. 


mixture  technique 
characteristic  functions 
geometric  transform 


A&G1.  ANCKER,  C.J.,  Jr.,  and  QATAR  IAN,  A.V.,  "The  Distribution  of  Rounds 
Fired  in  Stochastic  Duels,"  Naval  Research  Logistics  Quarterly, 
Vol.  11,  No.  4,  Dec.,  1964,  pp.  503-327. 

(Also,  Systems  Development  Corporation,  Santa  Monica,  CA. . Document 
SP-1017/004/00,  4 March,  1964,  35pp.,  DDC  No.  AD-433  764.) 

FM-CRIFT)-  Fandom  ammunition  supply  (contains  fixed 
. fxft  ) supply  case ) 

General  Solutions : P(h),  P(h),  P(N«n|H) 

E[n|h),  BtN2  | Hi, 

PlKlJi^lH),  P[N-n|H],  P[N » n] 

Examples ; Distribution  of  N 
(no. of  rounds  fired) 

(1 ) Geometric 

(2 ) Infinite  supply 

(3)  Finite,  fixed  supply 

FD-  (l)  CRIFT  - (a)  random  ammunition  supply 
(b)  fixed  ammunition  supply 

General  Solutions : F[NA*n|A],  P[Na>e0|a],  P[Na ■ n | AB] , 

P[NA-n|B],  P[»A  - n) 
e[na|a],  e[i£|a] 

Marginal  increase  in  P[A]  if  ammunition  supply  is  Increased. 
Examples : 


A 

B 

pmf  of  Na 

IFT  pdf 

pmf  of  N^ 

EFT  pdf 

(1) 

Geometric 

ned 

Geometric 

ned 

(2) 

Binomial 

ned 

Infinite 
(no  limitation) 

ned 

(3) 

Fixed  (constant) 

ned 

Fixed  (constant) 

ned 

FD  - (2)  FD'T  - (a)  random  amunitlon  supply 
(b)  fixed  ammunition  supply 
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(1)  Geometric 

(2 ) Infinite  supply 

(3)  Fixed  (constant) 


Geometric 
Infinite  supply 
Fixed  (constant) 


mixture  technique 
characteristic  functions 
number  theory 


ANCKER,  C.J.,  Jr.,  and  GAFARIAK,  A.V.,  "The  Distribution  of  the  Time- 
Duration  of  Stochastic  Duels,"  Naval  Research  Logistics  Quarterly, 
Vol.  12,  Nos.  3 fi*  4,  Sept. -Dec.,  1965*  pp.  275-294. 

(Also,  System  Development  Corporation,  Santa  Monica,  CA.,  Document 
SF-1017/OOT/OO,  10  August,  1964,  29  pp.,  DDC  No.  AD-606  169.) 

PM  - CRIFT  ) (l)  random  time  limitation 
• FIFT  j (2 ) fixed  time  limitation 

General  Solutions : P[H],  pdf  time  to  hit,  pdf  time  for  no 

hits,  all  moments  of  tvo  preceding  pdf's 
(special  case,  no  time  limit),  pdf  of  total 
time  to  completion 


Example : 


IFT  pdf  pdf  - time  limit 

Erlang  (2 ) ned 


I 


A 4»  G2 . (cont'd) 

FD-  (l)  CRIFT 

- (2  ) FIFT 


U) 

<b) 


(a) 

(b) 


random  time  limitation 
fixed  time  limitation 

random  time  limitation 
fixed  time  limitation 


General  Solutions : pdf  of 


taIa, 


pdf  of  T.-  | AB,  all 


AB 

moments  of  the  preceding  pdf's  (special 
case,  no  time  limit),  total  time  pdf. 


Example : 


A's  IFT  pdf 

B's  IFT  pdf 

pdf  - time  limit 

ned 

ned 

ned 

Erlan  g (2 ) 

Erlang  (2 ) 

ned 

ned 

ned 

fixed 

Erlang  (2) 

Erlang  (2) 

fixed 

fixed 

(integer  X B's 

IFT) 

fixed 

ned 

fixed 

(integer  x B's 

IFT) 

fixed 

fixed 

fixed 

fixed 

none 

(infinite 

(1) 

(2) 

(3) 

00 

(5) 

(6) 

(7) 


mixture  technique 
characteristic  functions 
number  theory 
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A&W1.  ARCKEK,  C.J.,  Jr.,  and  WILLIAMS,  Trevor,  "Some  Discrete  Processes 
in  the  Theory  of  Stochastic  Duels,"  Operations  Research.  Vol.  13, 
Wo.  2,  Mar. -Apr.,  1965,  PP*  202-216. 

(Also,  System  Development  Corporation,  Santa  Monica,  CA.,  Document 
SB-101T/002/00,  13  August,  1963,  28  pp.,  DDC  Ko.  AD-t«20  51*0 

FD  - (l)  FIFT 

General  Solutions : P[A],  P[B],  P[AB] 

Examples : 

(1)  A's  IFT  an  integral  multiple  of  B's  IFT 

(2)  B's  IFT  an  integral  multiple  of  A's  IFT 

Solution  curves 

FD  - (2 ) Equal  FIFT  - probability  of  a neer  miss  is  included.  A near 

miss  causes  a displacement  and  the  loss  of  one 
firing  turn 

General  Solutions : P[A],  P[B]j  P[AB] 


Solution  curves 


mixture  procedure 
number  theory 

stochastic  difference  equations  for  (2 ) 
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Be  1.  BARFOOT,  C.B. , "The  Lanchester  Attrition  Rate  Coefficient:  Some 
Comments  on  Seth  Bonder's  Paper  and  a Suggested  Alternate  Method," 
Operations  Research,  Vol.  17,  No.  5,  Sept. -Oct. , 1969. 

# 

FM-  Markov  dependent  states  >2  - IFT  constant  hut  dependent  on 

current  state 

Kill  probabilities  dependent  on  current  state 
Initial  conditions  may  be  varied 

General  Solution:  E[T]:  in  matrix  form 

Example:  Numerical 


Markov  chain  theory 
matrix  algebra 


l 

/ 

< 

* 

1 


l 


Ba2.  BARFOOT,  C.B.,  "Markov  Duels,"  Operations  Research,  Vol.  22, 

No.  2,  Mar. -Apr.,  1974,  pp.  318-330. 

(Also,  "Stochastic  Duels  in  Which  Each  Contestant's  Shots  Form  a 
Markov  Chain,"  OR- 69.  5-th  International  Conference  on  Opera- 
tions  Research,  Venice,  Italy,  23-27  June,  1969,  ed.  by  John 
Lavrence,  Tavistock  Publ. , London,  1970,  pp.  223-234.) 

(Also,  Master's  Degree  Thesis  submitted  to  the  Department  of 
Operations  Research,  George  Washington  University,  Washington,  D.C.) 

(Also,  "Stochastic  Duels  With  Markov  Dependent  Kill  Probabilities," 
Center  for  Naval  Analyses  Working  Paper,  Arlington,  VA.,  undated, 

48  pp. ) 

FM  - FIFT  - Markov  dependent  states  > 2 

General  Solutions : pmf  (number  of  rounds  to  a kill  or  time  to 

kill);  in  matrix  form  E(n),  V(n) 

FD  - FIFT  - Markov  dependent  states  > 2 

General  Solutions:  P[A],  P[B],  P[AB]  in  matrix  form 

Example : Two-state  case  FM-FI7T  result 

FD  - FIFT  - M kov  dependent  states  >2.  A fires  y rounds  first 
(random  initial  surprise),  pmf (Y )-  geometric 

General  Solutions : P[A],  P[B],  P[AB]  in  matrix  form 


Numerical  example 
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Ba  2.  (cont'd) 

Markov  chain  theory 
matrix  algebra 


Ba  3.  BARFOOT,  C.B.,  "Some  Anti-Armor  Models  Used  in  U.  S.  Marine  Corps 

• Harming  Studies,"  NCQAG  CNA,  Arlington,  VA.,  NATO  Conference, 

26-50  A-igust,  1971*,  Hmich,  Geraany,  19  pp. 

FD  - FIFT  - Markov  dependent  states  > 2 

A fires  busts  of  constant  length,  with  a constant  time 

• between  rounds  and  a constant  time  between  bursts 

B fires  with  no  burst,  just  constant  time  between 
rounds 

General  Solutions : PfA],  P[B),  P[AB]  in  matrix  form;  asymptotic 

approximations  to  general  solutions 

Same  as  above  - but  A fires  y rounds  first  (surprise) 

(c.)  y a constant 

(b)  T a geometric  random  variable 

General  Solutions : P[A],  PlB],  P[A3]  in  matrix  form;  asymptotic 

approximations  to  general  solution 

FM  «■  CR3FT  - Markov  dependent  states  > 2 

IFT's  are  ned  and  dependent  on  current  state 
(different  parameter  for  each  state) 

General  Solution;  h(t) 

FD  - CRIFT  - Markov  dependent  states  > 2 

Each  duelist's  IFT's  are  ned  and  dependent  on  current 
state  (different  parameters  for  each  state) 

General  Solution:  F[A],  P[B] ; matrix  double  Integral  and  a closed 

solution  by  similarity  transformation 

FD  - CRIFT  - Markov  dependent  states  > 2 

Each  duelist's  IFT's  are  ned  and  state  dependent: 

(l)  A has  a fixed  time  to  fire  first  (surprise) 

£ (2 ) A has  a random  surprise  time  with  ned 

distribution 
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Bh  1.  BHASHTAM,  N.,  "Stochastic  Duel  With  Several  Types  of  Weapons,” 

Defence  Science  Journal  (India),  Vo!.  17,  Ko.  2,  April,  1967, 
pp.  115-118.' 

(Also,  Defence  Science  Laboratory  Report,  Delhi-6,  India,  9 pp. ) 

(Also,  Bh  3,  pp.  67-72 . ) 

FM  - CRIFT  - several  weapons  firing  simultaneously  and  independently, 
each  with  individual  ned  IFT's  and  kill  probabilities. 

General  Solution:  h(t) 

FD  - CRIFT  - each  side  with  several  (different  for  each  side)  weapons, 
firing  simultaneously  and  independently.  Each  weapon  with 
individual  ned  IFT's  and  kill  probabilities. 

General  Solution:  P[A3,  F[B] 

Mean  and  variance  of  number  of  rounds  of  each  type 
weapon  fired  to  a kill. 

differential  difference  equation  technique 
z transforms 

Bh  3 uses  elementary  methods  (much  simpler  but  cannct  get  last  results 
above. 


Bh  2.  BHASHTAM,  N.,  "Stochastic  Duels  With  Pattern  Firing,"  Advancing 

Frontiers  In  Operational  Research,  Proceedings  of  the  International 
Seminal’  on  Operational  Research,  New  Delhi,  India,  7-10  August,  1967, 
Ed.  by  H.S.  Subba  Rao,  N.K.  Jaiswal,  and  A.  Ghosal,  Hindustani 
Publishing  Corporation  (India),  1969 , pp.  151-164-. 

(Also,  Bh  3,  pp.  100-114) 

FM  - CRIFT  - two  weapons  fired  alternately,  each  firing  a fixed  number 
of  rounds  (different  for  each  weapon)  with  a different 
IFT  pdf  and  a different  hit  probability. 

General  Solution:  LT  h(t) 

FD  - CRIFT  - each  side,  two  weapons  fired  alternately,  each  firing  a 
fixed  number  of  rounds  (different  for  each  weapon)  with 
a different  IFT  pdf  and  a different  hit  probability 

General  Solution:  LT  time  to  a kill  by  A 

P[A],  P[B] 
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(cont'd) 
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Examples : 

(a)  Each  side  fires  one  round  each  with  two  weapons  and 
ned  IFT's 

(b)  A has  one  weapon,  one  round  each 
B has  two  weapons,  one  round  each 

(c)  A has  two  weapons,  one  round  each 
B has  one  weapon,  one  round  each 

Numerical  example,  curves 

Bh  3 - different  numerical  examples  and  curves 


supplementary  variable  technique 
differential  difference  equations 
geometric  transforms 
Laplace  transforms 


Bh  3-  BHASHIAM,  N.,  "Stochastic  Duels,"  Hi.D.  Thesis,  University  of  Delhi, 
Delhi,  India,  May,  1969,  186  pp. 

Pour  sections  in  this  thesis  have  not  been  published  in  the  open 
literature.  Only  these  are  annotated  here.  The  other  sections 
are  cross  referenced  in  the  appropriate  document. 

(l)  "Stochastic  Duels  With  Only  Pooled  Interfiring  Time  Distribution 
Known,"  pp.  75-8U. 

FD  - CRIST  - Probability  that  either  fires  next,  given  either  fired 
last,  is  probablistic . Therefore,  the  sequence  of 
firings  is  Markov  and  independent  of  times  between 
rounds.  IFT's  are  different,  depending  on  who  fired 
last. 

General  Solution:  LT  (l)  time  to  win  by  A (i.e.,  during  duel) 

(2)  time  to  win  by  B (i.c.,  during  duel) 

Numerical  example  - curves 

differential  difference  equations 
supplementary  variables 
geometric  transform 
Laplace  transforms 

(2)  "Stochastic  Duels  With  Burst  Fire,"  pp.  87-100 


I 

I 


Bh  3. 


(cont'd) 


FM  - CRIFT  - Bursts  of  fixed  length  fired  with  CREPT  and  with  contin- 
uous FV  between  bursts 

General  Solution:  LT  h(t) 

FD  - CH3FT  - Bursts  of  fixed  length  fired  with  CREFT  and  with  con- 
tinuous RV  between  bursts 

General  Solution:  LT  P[A),  P[B] 

Examples : (l)  all  times  ned,  burst  length  very  large 

(2 ) all  times  ned,  A large  burst  size,  B single 

round  fire  (no  bursts  J 

Numerical  results  - curves 


differential  difference  equations 
supplementary  variable  technique 
geometric  transform 
Laplace  transforms 

(3)  "Stochastic  Duels  of  Limited  Time  Duration  and  Finite  Ammunition 
Supply,"  pp.  1^1-153. 

FD  - CRIFT  (ned  both  sides ) - Both  fixed  ammunition  limitation  - 
time  limitation  a continuous  RV  with  ned  pdf.  A 
draw  occurs  if  time  runs  out  or  both  run  out  of 
ammunition 


General  Solution: 


Examples : jlj 

(1) 

(M 


P[A],  F[B] , P[AB] 
unlimited  time 

B unlimited  ammunition  - solution  curves 
B unlimited  ammunition  - unlimited  time 
both  unlimited  ammunition 


differential  difference  equations 
Laplace  transforms 

(4)  "Stochastic  Duels  with  Repairable  Weapons,"  pp.  153-169. 

FM- CRIFT  - (a)  Fixed  limited  ammunition 

(b)  ammunition  limitation  a discrete  RV 


Time  to  failure  of  weapon  pdf  is  ned.  Repair  time 
is  CFV 

General  Solution:  LT  h(t)  and  time  to  run  out  of  ammunition 

FD- CREPT  - (a)  Fixed  limited  ammunition  (both) 

(b)  Ammunition  limitation  a discrete  RV 


B-21 


Bh  3-  (cont'd) 


t 


I 


Time  to  failure  of  weapons  are  pdf's  with  ned's. 
Repair  times  are  CRV'b.  Contestant  under  fire  during 
repair  time. 

General  Solutions : P[A],  P[B],  P[AB] 

Particular  Cases : 

PM  !*  FD  - (l)  A unlimited  amrmmition ; B fixed 
ammunition  limit 

FM  4»  FD  - (2 ) Both  unlimited  ammunition 

Examples : (l)  Unlimited  ammunition,  IFT's  ned 
Repair  times  ned 

(2)  Unlimited  ammunition,  IFT's  ned 

A repair  time  ned;  B failure  free 

Solution  curves 

differential  difference  equations 
supplementary  variable  technique 
geometric  transforms 
Laplace  transforms 


Bh  4.  BHASHYAM,  N.,  "Stochastic  Duels  With  Single  Shot  Kill  Probability 
Varying  As  A Function  of  Inter-Firing  Time  Interval,"  Defence 
Science  Laboratory,  Delhi-6,  India.  Draft  - Private  Communication, 
Spring  1970,  10  pp. 

FM-CRIFT  - Fixed  ammunition  limitation 

- Kill  probability  a function  of  IFT 

General  Solution:  LT  h(t) 

FD  - CRIFT  - Fixed  ammunition  supplies,  both  sides 

- Kill  probabilities,  functions  of  IFT's 

General  Solutions:  Integrals  of  LT  P[A],  P[B],  P[AB] 

Example : Infinite  ammunition  supplies 
IFT ' 8 are  ned 

Kill  probabilities  a negative  exponent ial  function 
of  IFT's 

differential  difference  equations 
geometric  transforms 
Laplace  trails  forms 
supplementary  variable  technique 
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Bh  5.  BHASHTAM,  N.,  "Stochastic  Duels  With  Bound  Dependent  Kill  Probability 
and  General  Inter- Firing  Times/'  Defence  Science  laboratory,  Delhi- 
6,  India.  Draft  - Private  Communication,  Spring  1970,  15  PP* 

FM-  CHUT  - Ammunition  limitation:  (a)  fixed,  and  (b'  random 
Round-dependent  hit  probabilities 

General  Solutions:  LT  h(t) 

FD-  CRIFT  - Ammunition  limitation:  (a)  fixed,  and  (b)  random 
Round-dependent  hit  probabilities 

General  Solutions : P[A],  P[B],  P[AB];  also, 

(1)  A fixed  ammunition  limit 
B infinite  supply 

(2 ) Both  have  infinite  supply 

Examples:  (l)  both  fixed  annunitian  supply,  ned  IFT’s  % 

(solution  in  terms  of  unspecified  hit  probabilities ; 

(2 ) both  fixed  ammunition  supply  and  both  general 

Erlang  IFT's  (solution  in  terms  of  unspecified  hit 
probabilities ) 

differential  difference  equations 
special  discrete  transforms 
geometric  transforms 
Laplace  transforms 
supplementary  variable  technique 


Bh  6.  BHASHYAM,  N.,  ''Stochastic  Duels  With  Non-Repairable  Weapons,"  Naval 

Research  Logistics  Quarterly,  Vol.  17,  No.  1,  March,  1970,  pp.  121-129- 

(Also,  Defence  Science  Laboratory  Report,  Delhi-6,  India,  undated, 

13  PP* ) 

(Also,  Bh  3,  PP*  169-181.) 

FM  - CRIFT  (ned  only)  - Limited  ammunition,  failure  wane  weapons  with 

u limited  replacement  stock  (failure  times  are 
ned) 

General  Solutions:  LT  h(t) 

LT  time -to- failure  (weapons  or  ammunition  supply) 
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Bh  6.  (ccnt'd) 


FD  - CRIFT  (ned  only)  - Limited  ammunition,  failure  prone  weapons  with 

limited  replacement  stock  (failure  times  are 
ned) 

General  Solutions : P[A],  P[B],  P[AB] 

Also,  same  with  unlimited  ammunition  and  replace* 
ment  stock  a discrete  RV 

Examples : (l)  both  unlimited  ammunition,  A limited  fixed  weapon 
supply 

numerical  illustration  with  curves 

(2 ) from  Bh  3 - both  unlimited  ammunition,  both 
geometrically  distributed  number  of  Weapons 

differential  difference  technique 
Laplace  transforms 


Bh  7. 


BHASHIAM,  N.,  "Stochastic  Duels  With  Lethal  Dose,”  Naval  Research 
Logistics  Quarterly,  Vol.  17,  No.  3,  Sept.,  1970,  pp.  397-1*05.  ~ 

(Also,  Bh  3,  pp.  11^-125.) 

FM-  CRIFT  - Multiple  hits  to  a kill  (P.) 

(a ) fixed  R 

(b)  R a discrete  RV 

General  Solution:  LT  h(t) 

FD-  CRIFT  - Multiple  hits  to  a kill  (r) 

(a)  fixed  R 

(b)  R a discrete  RV 

General  Solutions : P[A],  P[b3 

Examples : (l)  TFT  ned  for  both,  R fixed  for  both 

Solution  curve* 

(2 ) IFT  ned  for  both,  R a geometric  RV  for  both 

Solution  curves 

geometric  transforms 
differential  difference  equations 
Laplace  transforms 

Bh  3 gives  a much  simpler  derivation 
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Bh  8.  BHASHYAM,  N. , "Stochastic  Duels  With  Correlated  Fire,"  Metrlka, 

Vol.  20,  Ho.  i,  February,  1973,  FP*  17-24. 

(Also,  Bh  3,  pp.  125-137.) 

FM-  CRIFT  - Two  weapons  with  different  HT' s and  different  hit 

probabilities.  The  probability  of  firing  a given  weapon 
on  the  next  round,  given  a particular  weapon  was  fired  on 
the  last  round,  is  fixed.  This  leads  to  a correlation 
between  sequences  of  weapons  fired. 

General  Solution:  LT  h(t) 

FD  - CRIFT  - Two  weapons  with  different  IFT * 8 and  different  hit 

probabilities.  The  probability  of  firing  a given  weapon 
on  the  next  round,  given  a particular  weapon  was  fired 
on  the  last  round,  is  fixed.  This  leads  to  a correlation 
between  sequences  of  weapons  fired. 

General  Solution:  P[A],  P[B] 

Examples:  ned  IFT's  for  both 

differential  difference  equations 
geometric  transforms 
Laplace  transforms 
supplementary  variables 
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Bh&Si  1.  BHASHXAM,  N.,  and  SINGH,  N.,  "Stochastic  Duels  With  Varying  Single 
Shot  Kill  Probabilities,"  Operations  Research,  Vol.  15,  No.  2, 

Mar. -April,  1967,  pp.  2 53-2 4^ . ”” 

(Also,  Defence  Science  Laboratory  Report,  Delhi-6,  India,  November, 
1966,  19  pp. ) 

FD-CRIFT  (ned  only) -Fixed  ammunition  limitation 

- Kill  probability  is  a function  of  round 
number 

General  Solution:  (l)  P[A],  P[B],  P[AB] 

(2 ) B has  infinite  ammunition  supply  P[A] , 

P[B] 

(3)  A and  B have  infinite  ammunition  supplies 
P[A],  ?[B]  developed  separately  without  LT 

Examples ; (l)  fixed  kill  probabilltes ; B infinite  ammunition 
supply;  also,  both  infinite  ammunition  supply 


(2)  P(n>  - 


(3)  p(n)  » (l-Q^1); 

oc,S  parameters 

differential  difference  equations 
special  discrete  transforms 
Laplace  transforms 


p*(m) 


1 

m t l 


p' (m)  « (l  - Pm) 
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Bo  1.  BONDER,  Seth,  "The  Lenchester  Attrition-Rate  Coefficient,"  Operations 
Research,  Vol.  15,  No.  2,  March-April,  1967,  pp.  221-232. 

FM  - PITT  - Markov  dependent  hit  probabilities  (dependent  on  states ) 
three  states 

Fixed  number  of  multiple  hits  required  to  kill 

General  Solution;  pof  of  N - number  of  rounds  fired  to  give  a 

fixed  number  of  hits 

E[N] 

combinatorial  arguments 


Bo  2.  BONDER,  Seth,  "The  Mean  Lanchester  Attrition  Rate,"  Operations 
Research,  Vol.  18,  No.  1,  Jan. -Feb.,  1970,  pp.  179-181. 

General  Solution:  If  T is  RV,  time  to  a kill  in  situation 

of  Bo  1,  gives  E[T] 

elementary  probability  arguments 


M 
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Fi  1.  FINLEY,  David  R. , "A  Theoretical  Study  of  Round- to-Round  Correlation 
In  Gunnery,"  Cornell  Aeronautical  Laboratory,  Inc.,  Buffalo,  New 
York,  Internal  Research  Report,  WA- 86-184,  Nov.,  1988,  33  pp. 

(Also,  Master  of  Arts  Thesis  at  the  American  University,  10  May,  1968, 

33  PP>  Available  at  University  Microfilms,  Inc.,  Ann  Arbor,  Michigan, 

No.  ML647) 

FM-  CRIFT 

I miss 

hit,  not  killed 
hit,  killed. 

P[hit]  ■ p,  a constant;  P[miss]  » 1-p;  and  P[kill|hit]  ■ 0,  a 
constant.  However,  p may  depend  on  any  or  all  of  the 
previous  rounds  fired.  General  results  are  derived  which  dc 
not  apply  directly  to  FM  or  FD  because  more  than  one 
killing  round  is  allowed.  However,  these  results  can  be 
adapted  to  FM  or  FD.  Verv  generally,  positive  correlation 
is  defined  as: 

P[hit  on  i-th  rcxind  | hits  on  specified  previous  rounds] 

> P[hit  on  i-th  round  | miss  on  at  least  one  of  the 
~ specified  previous  rounds  and  miss  on  all  other  rounds]. 

Note:  (a)  any  hit  may  be  a kill 

(b)  this  is  not  the  same  as  the  usual  definition  of  cor- 
relation and  the  results  do  not  apply  to  ordinary 
correlation. 

General  Solution:  Positive  correlation  decreases  kill  probability, 

compared  to  the  case  where  all  trials  are  inde- 
pendent. For  usual  definition  of  correlation, 
positive  correlation  does  not  necessarily  decrease 
kill  probability. 

no  overkills  (i.e.,  first  kill  terminates  process) 
dependence  iB  Markov 

positively  correlated  if  correlation  coefficient  > 0; 
distribution  of  N,  E[N],  V[n] 

Example : cf  of  h(t)  for  general  IFT  and  E[T],  V[T] ; 
also  ned 

FD  - A is  PM  as  above,  B is  FM 
Both  IFT*  8 ned 

General  Solution : P[ A] 


Example : 

Results : 


» 


set  theory 
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Ft  1.  FRIEDMAN,  Yoram,  "A  Model  for  the  Determination  of  Optimal  Inter- 
firing  Times,"  (unpublished  ms.),  Faculty  of  Management,  Tel  Aviv 
University,  Tel  Aviv,  Israel,  July,  1976,  8 pp. 

FM  - Fixed  time  limit 

• Kill  probability  a monotone  non-decreasing  function  of  TFT 
- IFT  is  continuous  but  greater  than  a certain  minimum  and  at 
choice  of  f irer 

For  a fixed  number  of  rounds  fired  (n)  [ < the  maximum 

possible  in  time  limit]  proves  that  optimal  kill  policy  is 
to  continually  decrease  the  IFT's,  but  use  up  all  the  time 
available . 

Solution : Gives  an  algorithm  for  finding  the  t‘s.  Also  shows 
how  to  find  the  best  n. 


Lagrange  multipliers 


Ga&Sl.  GARG,  R.C.,  and  SINGH,  N.,  "A  Stochastic  Duel  in  a Hunter-Killer 
Game  - I,"  The  Symposium  on  Operations  Research,  No.  42,  1970, 
pp.  185-192 


FD  - CRIFT  (both  ned) 


5 states 


(1)  combat  - time  is  ned  v 

(2)  no  contact  - time  is  ned  ' all  CRV's 

(5)  reclose  & continue  time  is  J 
general  pdf 


General  Solutions : LT’s  of  time-functions  of  various  states;  in 
“ particular  that  A or  B has  von  at  time  t 

after  each  has  expended  a certain  number  of 
rounds 

Examples : (l)  reclose  and  continue  time  is  ned 
(complete  solution) 

(2)  same  as  (l),  but  FM  only 


differential  difference  eauations 
Laplace  transforms 
supplementary  variable  technique 
geometric  transforms 


( 

L 

r 
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Gr  1.  GROVES,  Arthur  0.,  "The  Mathematical  Analyst*  of  a Simple  Duel," 
Ballistic  Research  Laboratory  Report  Ho.  1261,  Aberdeen  Proving 
Grounds,  Md. , August,  1964,  23  pp. , DEC  No.  AD-609  195* 

PD  - 7IFT  - A starts  with  a fixed  time  advantage  (time  he  fires 
before  B does ) 


t 

t 

f I 


General  Solution;  in  matrix  form 

(1)  state  probabilities  after  m cycles 

(2)  P[A],  P[B],  P[AB] 


Numerical  examples 

L 


f 

I-  Markov  chains 

[ number  theory  - matrix  manipulations 

• geometric  transforms 


£ 


ii 

L. 
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HARRIS,  Terrell  J. , Many  Versus  Many  Stochastic  Duels,"  Cayvood 
Schiller  Associates  Report,  Chicago,  Illinois,  Fall  1967,  18  pp. 

FD  - CRIFT  - A general  IFT 
B ned  IFT 

General  Solutions:  for  P[Al*  P[B]  in  LT  and  geometric 

transforms  for  mean  number  of  rounds  fired 
to  a victory  for  A 

Examples  : (l)  A's  IFT  > Erlang  (k) 

(2)  A's  IFT -mixture  of  ned's 


renewal  theory 
geometric  transforms 
Laplace  transforms 


FM  - CRIFT-  fixed  ammunition  supply  at  start 

- ammunition  replenishment  of  a fixed  amount  at  times 
with  pdf  (ned ) 

- TCF  is  CRV 


l 

r 

£ 

jr 

l 

\ 

i 

l 

f 

\ 
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General  Solution:  LT  of  h(t) 


Special  cases : 


(1)  no  replenishment,  LT  of  h(t) 

(2)  no  replenishment , initial  supply  a discrete 
RV  (in  BH  3 only),  LT  of  h(t) 

(3)  no  replenishment,  unlimited  ammunition, 

LT  of  h(t) 

(4)  sero  flight-time,  LT  of  h(t) 

(5)  sero  flight-time,  no  replenishment  and 
initial  supply  is  fixed,  LT  of  h(t) 

(6)  sero  flight-time,  no  replenishment  and 
initial  supply  is  a discrete  RV,  LT  of 

(7)  flight-time  sero,  no  replenishment,  un- 
limited ammunition,  LT  of  h(t) 


FD  - CRIFT  - fixed  ammunition  supply  at  start 

- ammunition  replenishment  of  a fixed  amount  at  times  vith 
pdf  (ned) 


V 

- TOF  is 

CRV 

General  Solutions 

P[Aj,  P[Bj,  P(AB) 

Examples  : ) 

[TOT'S  ned  " unlimited  ammunition 

t 

I 

y 

(2)1 

I A^s  TOT  ned  * eomunition,  sero  flight 

time  for  B 

curves  in  BH3 

ft 

(3) 

IFT's  ned -both  sero  flight-time,  unlimited 
ammunition 

w 

00 

IFT's  ned -both  flight-times  sero  (in  Bh3) 

- B unlimited  ammunition 

differential  difference  equations 
geometric  transforms 
Laplace  transforms 
supplementary  variables 


Ki  1.  KIMBLETGN,  Stephen  R. , "Attrition  Rates  for  Weapons  With  Harkov- 
Dependent  Fire,"  Operi.y \ons  Research.  Vol.  19,  No.  3,  May-June, 

1971,  pp.  698-706. 

PM  - FIFT  - Markov  dependent  hit  probabilities  (dependent  on  states ) 

3 states,  random  -.lumber  of  multiple  hits  to  kill 

General  Solutions : Laplace  transform  of  N (number  of  rounds  to 

a kill),  E[N),  V[N],  pmf  of  N 

If  T is  time -to- kill,  Laplace  transform  of 
pmf  of  T,  E[T),  V[T] 

H 

Markov  chain  theory 

renewal  theory  1 

difference  equations  i 

Laplace  transforms  j 

geometric  transforms  j] 
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KW&B1.  KWCN,  T.Y.,  and  BAI,  D.S.,  "Stochastic  Duels  With  Multiple  Hits 
and  Limited  Ammunition  Supply,"  Korea  Advanced  Institute  of 
Science  Report,  Seoul,  Korea,  April,  1978,  23  pp. 


» 


FM-  GRIFT  (l)  Fixed  nunber  of  hits  to  a kill 

(2  ) Random  number  of  hits  to  a kill  - geometrically 
distributed,  i.e.:  Parameter  - Pr [kill  | hit] 

(3)  Fixed  nunber  of  hits  to  a kill;  ammunition  limit- 
ation a discrete  RV 


(k)  Random  nunber  of  hits  to  a kill  (geometrically 
distributed);  ammunition  limitation  (geometric- 
ally distributed) 

(3)  Fixed  number  of  rounds  fired  simultaneously 

(pattern);  fixed  probability  of  a pattern  hitting, 
each  round  in  pattern  which  hits  has  fixed 
probability  of  a kill 

(6)  Limited  ammunition  a discrete  RV;  fixed  nunber  of 
rounds  fired  simultaneously  with  fixed  probability 
of  a pattern  hitting  and  each  round  in  a pattern 
kills  with  a fixed  probability 


General  Solutions  : LT  h (t ) 


FD  - CRIFT  - Both  sides  same  as  FM  (l),  (2)  and  (3) 

General  Solutions : P[A],  P[B],  P[AB] 

Examples : (a)  fixed  number  of  hits  to  a kill;  ammunition  limit 
a geometric  discrete  RV;  IFT’s  ned 

(b)  geometric  number  of  hits  to  a kill;  geometric 
distribution  for  nunber  of  rounds,  IFT's  ned 

- curves 

(c)  pattern  firing  (3)  above,  with  IFT's  ned 

(d)  pattern  firing  (6)  above,  with  geometric  distri- 
bution for  number  of  rounds,  IFT's  ned 

- curves 

mixture  technique 
Laplace  transforms 

Case  (2)  also  derived  using  renewal  theory 


t 


t 
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M&A1.  MOHAN,  C.,  and  ARORA,  S.D.,  “On  a Problem  in  Naval  Defense," 
Operations  Research  , Vol.  12,  No.  2,  March  - April,  1963, 
pp.  19^-198. 

FM-CRI7T  (ned) 


2 states 


(l)  combat  - tine  length  is  ned 

(2  ) between  engagements  - time  length  general 


CRV's 


Although  ITT  la  ned,  only  firing  times  when  hits 
occur  are  considered,  so  these  events  are  ned  with 
a parameter  pr  (not  r) 


General  Solution;  LT  of  time  to  n hits 

Example : between  » betveen-engsgements  time  ned 
also,  expected  time  for  n hits 


differential  difference  equations 
supplementary  variable  technique 
Laplace  transforms 
geometric  transforms 
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N&J1.  XAGABHLJSHAXAM,  A. , and  JA3X,  O.C.,  "Stochastic  Duels  With  Damage," 
Operations  Research,  Vol.  20,  Xo.  2,  March- April,  1972, 
pp.  350-556. 

(Alto,  Defence  Science  laboratory  Report,  Delhi-6,  India,  undated, 
13  pp.) 


FM  - CRIFT  (l ) Amount  of  damage  per  round  is  a discrete  RV  with 
a jmf . Damage  is  Independent  round- to-round  and 
cumulative  until  total  is  a kill  (predetermined 
total  damage) 

(2)  Damage  per  found  is  a CRY  with  pdf;  otherwise, 
same  as  (l)  above 

(3)  Damage  is  time -homogeneous,  i.e. , damage  increase 
in  At  * p ( increase )At  + 0(At  ),  where 
Increase  is  a discrete  RV 


General  Solutions : cf  h(t) 

Example ; three  discrete  damage  states  for  (l),  IFT  ned 

FD  - CRIFT  - same  for  both  as  (l),  (2;  ) above,  plus 

(U)  Damage  states  are:  no  damage,  damage,  Kill- 
Given  damage,  the  amount  is  a CRY  with  a 
cumulative  upper  limit  ctv  iing  a kill. 

General  Solutions : P[Aj,  P[B] 

Examples:  (a)  three  discrete  damage  states  t'Ji  (l),  IFT’s 
ned 

(b)  for  (3)  above,  two  discrete  damage  ^ cates 


difference  equations 
geometric  transforms 
X renewal  theory 

differential  difference  equations 
characteristic  functions 


» 


t 


R* 81.  RUSTAGI,  J.S.,  and  SRIVASTAVA,  R.C.,  "Parameter  Estimation  in  a 
Markov- Dependent  Firing  Distribution,”  Operations  Research, 

Vol.  16,  No.  6,  Nov. -Dec.,  1968,  pp.  1222-1227. 

FM  - F2TT  - Two  Markov-dependent  states 

- Multiple  hits  (r)  to  get  a kill,  different  first 
round  hit  probability 

General  Solution:  pgf  (geometric  transform)  of  N (number 

of  rounds  to  get  r hits) 

MLE  of  the  three  parameters  involved 

Markov  chain  theory 
renewal  theory 
geometric  transforms 


St  1. 


SAV1B,  David,  "Asynchronous  Dodging  Duels,"  private  ccoaunication, 
19  August,  1970,  pp.  5^-69. 


FD  - FIFT  - equal  IFT's 

A • tarts  first,  B starts  later,  by  an  amount 

Each  can  either  hit,  miss  or  near-miss,  causing  1 
placement  and  the  loss  of  one  firing  time  (while 
remaining  vulnerable . ) 

General  Solutions : P[A] , F[B] 


Boolean  algebra 


< IFT 
1 dis- 


(3)  Equal  firing  times ; A fires  first  vlth  probability  p. 

A evacuates  after  k > 3 rounds,  evacuation  distance  RV 


General  Solution: 


P[A  Is  alive  time  T after  Initiating  evacuation 
alive  to  evacuate] 


A is 


Numerical  example 
combinatorial  methods 

Sc  2.  SCHODERBEK,  J. J. , "Some  Weapon  System  Survival  Probability  Models - 
II.  Random  Time  Between  Firings,"  Operations  Research,  Vol.  10, 

No.  2,  March- April,  1962,  pp.  168-179* 

FD  - CRDT  (both  ned) 

(l)  Kill  probabilities  and  rates  of  fire  are  continuous 
functions  of  time 


H 


numeric*!  Examples:  kill  probabilities  and  rates  of  fire  constant 
FD  - CRIFT  (both  ned) 

(2)  Same  as  (l),  but  A evacuates  at  time  t0  (A  la  vulnerable 
during  evacuation).  A moves  a distance  which  Is  Rayleigh 
distributed . 

General  Solution;  F[A  survivies  evacuation  | A Is  alive  to  start] 
and  P[A  survives  evacuation] 

differential  equations 


I Sr  a.  01. 


f 
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SRIVASTAVA,  S.S.,  and  GARG,  R.C.,  "Two-Sided  A/S  Warfare  With 
Limited  Ammunition  and  Range  of  Fire  Power,"  Traba.1  ob  Ae 
Bstadlstica  y de  Inveatleacion  Operative,  Vol.  XXIII,  Madrid, 
SpainTT^Tw^nS-I1!?^ 

FD  - CRIFT  - IFT'b  (both  ned)  - both  limited  ammunition 


it  states: 


(1) 

(2) 

(3) 

OO 


Seeking  contact  - time  is  general 

Closing  in -time  is  general  pdf 

Combat  - time  is  ned 

Release  from  combat  - time  is  ned 
(before  starting  (l)  again) 


All  CRV's 


General  Solutions : LT's  of  time  function  of  various  states  given 

above*  In  particular,  that  A wins  or  £ 
wins,  or  a draw  has  occurred  at  time  t, 
after  each  has  fired  a specified  number  of 
rounds 


Example : seeking  contact  and  dosing  in;  both  ned;  general 
solutions  as  given  above,  inverse  LT's  of  these 
functions  and  same  functions  at  t » • 


differential  difference  equations 
supplementary  variable  technique 
Laplace  transform 


trr^c 
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SRIVASTAVA,  S.S.,  GARG,  R.C.,  and  SINGH,  N.,  "A  Stochastic  Duel  Ir 
a Hunter- Killer  Game  - III,"  Cahiers  du  Centre  d* Etudes  de 
Recherche  Operatioonelle  (Belgium;,  Vol.  11,  No.  2,  July,  1969, 
pp.  104-111. 


FD  - CRIFT's  (both  ned) 


U states : 


(l)  Closing 


(2 ) Combat 

(3)  No  contact 

(b)  Reclose  S»  continue 


/ All  these  times  are  CRV’s  with  general  pdf’s 

) 


General  Solutions;  Bivariate  geometric  transform  of  LT  of 

time  functions  of  various  states.  In 
particular,  that  A (or  B)  has  won  at 
time  t after  each  has  expended  a specified 
number  of  rounds. 


differential  difference 
supplementary  variables 
geometric  transforms 
Laplace  transform 


i 


’ 

l 


equations 

technique 
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THOMPSON,  David  E.,  "Reliability  and  Mobility  in  the  Theory  of 
Stochastic  Duels, ' ftu*t  F,  Ch.  1,  pp.  575-612,  in  Development  of 
Analytical  Models  of  Battalion  Task  Force  Activities,  ed.  by  S. 
Bonder  and  R.  Farrell.  Report  No.  SRL  1957,  FR  70-1 (u),  Systems 
Research  Laboratory,  Dept,  of  Industrial  Engineering,  University  of 
Michigan,  September,  1970,  DDC  No.  AD 7146 77. 

(Also,  "Mobility  and  Reliability  in  the  Theory  of  Stochastic  Duels," 
Master's  Thesis,  Dept,  of  Industrial  Engineering,  University  of 
Michigan,  Ann  Arbor,  Michigan,  6 August,  1968.) 

(Also,  DSL  CE147  W?  68-4(u),  Defense  Systems  Laboratory,  University 
of  Michigan,  Ann  Arbor,  Michigan,  31  October,  1968. ) 

(Also,  "Reliability  and  Mobility  in  the  Theory  of  Stochastic  Duels," 
a chapter,  pp.  575-612,  in  Development  of  Models  for  Defence  Systems 
Hanning,  ed.  S.  Bonder  and  R.  Farrell,  Technical  Report  (u) 

SRL  2i Vf  TR  70-2,  Systems  Research  Laboratory,  Dept,  of  Industrial 
Engineering,  University  of  Michigan,  Ann  Arbor,  Micnlgan,  September, 
1970- ) 

FD  - CRIFT 

(l)  Weapon  life  time  a CRV 

(a)  no  withdrawal 

General  Solution:  P[A],  F[B],  P[AB]  in  terms  of  pdf's 

and  df 's 

Example;  IFT's  - ned,  weapon  lifetimes  ned 

(b)  withdrawal  at  failure  times 

General  Solution:  PU’J,  P[B],  P[AB]  in  terms  of  pdf's 

and  df 'b 

Example ; IFT's  • ned,  weapon  lifetimes  ned. 

Can  be  shown  to  be  same  as  A2(l)(a) 

(c)  withdrawal  at  next  firing  time  after  failure 

General  Solution;  F(A],  P[B],  PfAB]  in  terns  of  pdf's 

and  df 's 

Example ; IFT's  - ned,  weapon  lifetimes  - ned. 

Discussion  of  significance  of  failure  rates. 
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TX.  (cant'd) 


(2 ) Round-dependent  failures  - round  number  on  vhich  failure  occurs 
is  a discrete  RV 

(a)  no  withdrawal 

General  Solutions:  P[A],  P[B],  P[AB]  in  terns  of  pdf's 

and  df's 

Shown  to  be  sane  as  A1 

(b)  withdrawal  at  failure  tines  (no  delay); 
sane  as  A1  special  case 

(c)  withdrawal  at  next  round  after  failure  - fixed  kill 
probabilities 

General  Solutions:  P[Al,  P[B],  P[AB]  in  terms  of  pdf's 

and  df's 

Example?  IFT'e  ned 
FM  - CRIFT  (ned) 

(3)  Time-dependent  kill  probabilities 

General  Solution:  h(t)  in  terms  of  pdf's  and  df's 

FD-  CRIFT 's  (ned)  - Time-dependent  kill  probabilities 

General  Solutions ; P[Aj,  P[B],  P[AB]  in  terms  of  pdf’s 

and  df's 

Example ; 

p (t)  ■ £ — =—  , similarly,  for  B 

A <r8+vt2) 

(a,  r and  v,  constants ) 

8 

FM  - CRIFT  (ned) 

(It)  Time- dependent  km  probabilities  and  time- dependent  firing 
rates 

General  Solution:  h(t)  in  terms  of  pdf's 

and  df's 
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Tl.  (cont’d) 

P.(t)  r.(t) 

PD  - ••me,  end  special  case:  ■ k 

B B 

elementary  probability  arguments 
mixture  technique 

(3)  uses  stochastic  differential  equation  method 

T2.  THOMPSON,  David  E.,  "Stochastic  Duels  Involving  Reliability," 

Naval  Research  Logistics  Quarterly.  Vol.  .19,  No.  1,  March,  1972, 
pp.  145-148. 

Covers  some  but  not  all  of  Tl 

One  new  case,  2(c),  round-dependent  withdrawal  at  next  firing  time 
after  failure  and  probability  of  failure  on  any  given  round  is  a 
discrete  RV 

General  Solutions:  P[A],  P[B],  P[AB)  in  terms  of  pdf’s  and  df’s 


mixture  technique 

elementary  probability  arguments 


r 


Wl.  WILLIAMS,  Trevor,  "Stochastic  Duels  - II,"  System  Development  Corporation, 
I Santa  Monica,  CA. , Document,  SMOlT/OOjj/OO,  1?  Sept » , 1963,  6l  pp. 

DDC  Ho.  AD420-515- 

FM- CRIFT  - Cumulants  of  h(t ) (up  to  the  fourth)  derived  In  terms 
of  cumulants  of  3PT  pdf 

’ FD- CRIFT  - In  times-to-frlll  as  Erlang  (k)  functions 

General  Solution : Using  FM  above,  can  be  used  as  an  approximation 

to  any  FD  using  means  and  variances  of  IFT's 

Illustrations 

FD  - CRIFT  - Approximate  solution  in  terms  of  means  and  variances  of 
IFT's  and  kill  probabilities  (2  terms  of  an  Infinite 
series ) 

mixture  theory 

moment  generating  functions 

finite  calculus 


W2 . WILLIAMS,  Trevor,  "Stochastic  Duels  With  Homing,"  System  Development 
Corporation,  Santa  Monica,  CA.,  Document,  SB*1017/l06/00,  18  May, 

1965,  34  pp.,  DDC  No.  AD-617- 773. 

(Also,  "Stochastic  Duels  • III, " System  Development  Corporation,  Santa 
Monica,  CA.,  Document  SP-1017/006/0 0,  22  June,  1964,  72  pp. 

DDC  No.  AD-443- 754.) 

FD  - CRIFT  (both  ned  ) ^ 

Hit  probability,  (pn),  round-dependent  and  increasing,  with 

probability  of  a hit  on  n-th  round  (H  ) a discrete  RV  with 
pmf  a negative  binomial 

A 1 - , N 

p ■ hit  W ■ - where  n ■ E(Ny,  N the  RV,  round  on 
n 

which  hit  occurs 

This  causes  p^  < p < p„ 

General  Solution:  P[A],  P[B] 

- Curves  comparing  outcome  vs  pA,  p.  (instead  of  p^, , p^ ) 
for  various  parameter  values  A 

- MLE  of  p and  £ where  k is  a parameter  in  Erlang  (k) 


£ 


t 


mixture  theory 

moment  generating  functions 

geometric  transforms 
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W&A1.  WILLIAMS,  Trevor,  and  ANCKER,  C.J.,  Jr.,  "Stochaatic  Duels,"  Operations 

Research,  Vol.  11,  No.  5»  Sept  .-Oct.,  1963,  pp.  803-817. 

(Also,  System  Development  Corporation,  Santa  Monica,  CA.,  Document 

SP-1017/000/01,  20  March,  1963,  22  pp.,  DDC  No.  AD-l»00  637.) 

FD-  CRIFT 

General  Solution:  P[A],  P[B] 

Examples ; (l)  IFT's  both  ned  - curves 

(2)  IFT’s  both  Erlang  (2)  - curves 

FD  - CEIFT  - Both  fire  simultaneously  at  time  zero  ("classical  duel" ) 

General  Solution;  J?[A],  P[B],  F[AB] 

Example : IFT'j  both  ned  - curves 

FD  - CRIFT  - A fires  first  l/2  the  time,  B fires  first  l/2  the  time 
(tactical  equity) 

General  Solution:  P[A],  P[B) 

Example : IFT's  both  ned  - curves 

FD  - CRIFT  - For  a random  initial  period  of  time,  one  aide  or  the  other 
may  fire  vith  impunity  (i.e.,  the  other  aide  cannot  return 
fire)  - called  "Random  initial  surprise"  - time  advantage 
a CRV  - positive  values  are  A's  advantage,  negative 
values  are  B's  advantage 

General  Solution;  P[A],  P[B] 

2 

Example : IFT's  ned  and  surprise  time  N(o,c  ),  curves 


mixture  theory 
characteristic  functions 
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W&A2.  WILLIAMS,  G.  Trevor  (ed . by  ANCKER,  C.J.,  Jr.),  "Stochastic  Duels 
With  Displacements  (Suppression),"  U.S.  Army  TRADOC  Systems 
Analysis  Activity,  White  Sands  Missile  Range,  New  Mexico,  TRASAKA 
Technical  Memorandum  5-77,  March,  1978,  8 pp.,  DDC  No.  AO- 52116. 

FD-CRIFT  (both  ned) 

A miss  may  be  a complete  miss  or  a near  miss.  A near  miss  may 
either  be  a kill  or  cause  a movement  (suppression  of  fire)  which 
lasts  a time  which  is  a CRV  (both  ned; 

Contestant  is  vulnerable  during  a displacement  (suppression), 
but  cannot  return  fire 

General  Solution : P(A] 


difference  equations 
conditional  probabilities 


ZINGER,  A.,  # Concentrated  Firing  in  Many  Versus  Many  Duels," 
Universite  du  Quebec,  Montreal,  unpublished,  July,  1978,  27  pp 

FD  - Alternate  firing 

- Fixed  multiple  hits  to  a kill  required 

(different  for  each  side) 

General  Solutions : P(A  | B starts),  F(A  | A starts), 

P(B  | A starts ),  P(B  | B starts ) 


difference  equations 


